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We present a microscopic theory of superfluid *He, formulated using the 
overall experimental observations as input information. With the theory of a 
consistent basis, we answer all of the essential questions regarding He II. To 
conform to the overall key experimental observations, we show that in their 
coexistence at a higher T the superfluid and normal fluid atoms necessarily 
aggregate in separate large regions, which with falling T reduces to predomi- 
nantly superfluid aggregates. In a superfluid aggregate, termed a momenton, 
the helium atomic wave-packets primarily perform localized (harmonic) os- 
cillations about sites fixed on the momenton, facilitating excitation states of 
discrete energy levels. The predominant thermal excitation in the superfluid 
equilibrium is hence between these discrete levels, i.e. in the form of single 
phonons. Superposed to the coherent oscillations of the single atoms above, the 
momenton waves, comprising the in phase, simultaneous oscillations of all 
of the atoms in each momenton, give rise to the so-called second sound wave, 
which contribution to the total excitation and energy is however negligible. 
The superfluid atoms are strongly correlated, due to a many-quantum-atom 
correlation, having a binding energy —7.2 °K/atom as from experiments; from 
this we suitably deflne a superfluid bond, A(, ~7.2 °K/atom. At, is substan- 
tial, on the scale of total excitation energy (~ 15 °K) or of the binding energy 
of He I (~ 3 °K/atom), and inevitably dominates the phase stability of He II. 
Accordingly in the superfluid in quasi equilibrium between 0.6-2.17 °K, ther- 
mal excitation is predominantly due to the break-up of At, at an activation 
energy Af, + (Sv — 8.6 °K. This yields the experimentally observed A speciflc 
heat, and in turn the inelastic damped-Bragg scattering of neutrons at 
1.93 A~^. By embodying these concepts pertaining to static structure of the 
fluid, atomic bobnding and atomic dynamics, we establish equations of mo- 
tion for the superfluid atoms, both semi-classical and quantum mechanical, in 
a coordinate decomposition scheme. Based on their solutions, and the results 
from a corresponding statistical thermodynamic treatment, we have been able 
to quantitatively predict all the essential properties of He II, including the full 
excitation spectrum (dividing in g < 1.93 — (7/2, 1.93 — (t/2 < 1.93 and q = 1.93 
A~^), the two fluid fractions, the specific heat (dividing in the T^- and the 
\-Cv regions), the superfluid phase transition temperature Tx, and in a sep- 
arate paper II the critical velocity Vc, the total number of excitation states 
and the superfluid viscosity; all of which are in good overall good agreements 
with experiments. We have also evaluated the zero point energy, the poten- 
tial energy, the first order thermodynamic functions including U{T), A{T), 
S{T) and J(T) of which ^(r) can be satisfactorily compared with experimen- 
tal data; we also formally present the continuity and stability conditions of 
the superfluid, and the dynamics structure factor due to the A;, excitation. 
The microscopic scheme also facilitates a novel, consistent QCE superfluid- 
ity mechanism, presented in paper II. Also on the basis of the scheme, we 
show that the circulation qunataization is due to a circular atomic wave 
self-interference. As a whole, in terms of the theory proposed in this work 
not only are all the explanations regarding the essential properties of He II 
of a consistent theoretical grounding, but also are consistent with the concep- 
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tions that apply to common condensed matter systems and with the common 
laws of physics. The series of novel concepts evolving from this work can give 
significant impact also to the understanding of other superfluids. 



I. INTRODUCTION 

Liquid ^He exhibits properties which are altogether anonialous at temperatures below 
2.17 K, i.e. in the supcrfluid phase, He II. The two most striking ones arc its restraint from 
solidification down to zero K under ambient pressure^ and its superfluidity. Some 
of the anomalies, such as the large zero point energy^ and the volume expansion with falling 
T,^'^ can be traced immediately to the extreme chemical character of the ^He atom, which 
being light, small and inert. Many others, on the other hand, cannot be understood in this 
simple manner; the far-most of these are the superfluidity, the circulation quantization, ^^"^'^ 
the A specific heat about the phase transition, the anomalous excitation spectrum/^ and 
the anomalous thermal conductance^^ and melting curve. These properties can be seen to 
be largely dictated by the collective and correlation behaviour of the atoms in the condensed 
form and, their interpretations inevitably require firstly an adequate microscopic depiction 
of this condensed system. Since the discovery of the superfluidity of He II during 1911/'^ 
- 1938"'''*'* extensive investigations have been made on this system both experimentally 
and theoretically; although, experimentally direct probe of the excitations using neutron 
scattering technique, which informs a great deal about the microscopic dynamics, did not 
begin until the end of 1950's.^^ Many comprehensive reviews on this subject haven been 
given. 

Theoretically, the current understanding of He II is essentially based on the compilation of 
a few distinct models proposed by London,^'' Landau,^^ and Feynman^^ during 1938-1950's; 
a number of other studies^^"^^ have also historically contributed to the now accepted "sin- 
gle" model. With one thing in common, all these models assume that He II consists of two 
coexisting fluids, the normal fluid and the superfluid, and the two fluids interpenetrate on an 
atomic scale. The otherwise main concerns in each of the models have been the excitation 
scheme of the superfluid and a (corresponding) superfluidity mechanism. London'''*'^* sug- 
gested that of the tow fluids, the superfluid undergoes Bose-Einstein condensation (BEC), 
i.e. it consists of atoms that condense in the momentum space into a single degenerate 
state and in real space are motionless and non-interacting. Accordingly, the BEC implies a 
non-viscous and zero-entropy superfluid, and hence yields the siiperfluidity. This is despite 
the fact that, in his early view, London also foresaw the possibility that the quantum states 
of liquid He II would partly have the feature of the vibration states (phonons) of Debye 
in the theory of metals^^. As opposed to London, Landau^^ described He II as having 
two elementary excitations, phonons created mainly at a low T, and rotons at a higher T. 
Superfluidity will result if the velocity of a flow motion is lower than that of the particles 
from the elementary excitations. A landmark success of Landau's model is the excitation 
spectrum he phenomenologically designated in order to achieve a desirable (lambda) spe- 
cific heat; the key feature of the spectrum has been remarkably corroborated, firstly by the 
neutron scattering experiment of Larsson et al. in 1958.^^ However, Landau did not provide 
a microscopic scheme as to how the two elementary excitations occiir, and in particular, 
what actually a "roton" is. Furthermore, for the superfluidity Landau predicted a critical 
velocity which is at least two orders of magnitude too large than the experimental values, 
and contains no mechanism for channel width dependence. The two distinct models above 
were linked together by Bogoliubov^^ who suggested that the majority of superfluid is in 
BEC state in which, in contrast to London's BEC, the atoms interact weakly, and that 
the transition of the atoms from the BEC state to the quasi particle state produces the 
two types of elementary excitations of Landau. These two or three models are essentially 
phenomenological. Fcynman''^ in the 1950's, in a semi-microscopic description, suggsted the 
excitations (of both types) result from the superfluid, being the ground state, to the normal 
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fluid. Quantitatively based on variational principle for some trial wavefunction, where the 
microscopic picture for atomic motion and excitation is however essentially lost, he gave an 
improved prediction of the excitation spectrum (1954; Feynman and Cohen, 1956). At a mi- 
croscopic s(;alc;, Fcivnman briefly, qualitatively described the atoms as being confined in the 
cages built by neighboring atoms (- a picture held also in the present work), but primarily 
showed the "rotons" of Landau, representing the higher energy excitations, to be possible 
from a symmetry point of view. Taking Landau's critical velocity as a valid start, but to 
rescue it from the clearly large quantitative discrepancy, Feynman (1955) assumed there 
present even lower excitations, vortices; the creation of maximum vortices yields critical 
velocities giving the right magnitudes. Until 1990's, the widely adopted (excitation) model 
of He II is the up-dated compilation of the models by the four individuals, referred to below 
as the LLBF theory, is as follows:^^ the superfluid corresponds to a BE condensate, which 
becomes hundred percent at zero °K; the two elementary excitations result from the shut- 
tling of atoms between the superfluid - the BE condensate - and the normal fluid. Whilst 
the LLBF model assumes that, at zero °K, the BEC, being identified with the superfluid. 
is hundred percent, recent experiments have shown the BEC is only 9 ~ 10 percent. This 
conflict was circumvented by Clyde and Griffin in the 90's,^^ in proposing that the BEC 
fraction as well as the helium atom mass is modified due to an interaction of the superfiuid 
with the normal fluid atoms. This apparent improvement does not however remove the basic 
difficulties pertaining to the foundation of the LLBF model. Theoretical studies of liquid 
"^He in the past two decades have otherwise largely involved quantum simulations including 
variational MC^^ '"^^ and path integral MC,^^'^^ which essentially pertain to a more accurate 
representation of Feynman's (or, the LLBF) excitation model. More broadly, these studies 
are not yet apt to yield the microscopic information of the superfluid atoms, regarding con- 
flguration, motion and dynamics, and are even far remote from touching the microscopic 
scheme of superfluidity, and circulation quantization. The Green's function technique for 
treating many-body problems is not feasible today for treating such complex system as He 
II, in respect of its above ground propety in parcticular.^^'^^ 

As a whole, the LLBF theory encounters various difficulties which give rise to a variety of 
self-contradictions, and therefore does not represent an ultimately satisfactory theorization 
for the superfluid He 11. Its most significant difficulties may be outlined as follows (further 
discussions are pointed to each relevant sections). 1). The LLBF thermal excitation picture 
(when faithfully formally represented) inevitably predicts a broad peak in the dynamics 
structure factor, in direct contradiction to experimental observations (Sec. HA). 2). The 
non-viscous, non-interacting superfluid atoms, which although is sensibly not held in the 
dynamic model of the LLBF theory but is nevertheless a base for the LLBF explanation of 
superfluidity see item 4) below, directly contradicts with the fact that superfluid atoms 
are strongly correlated and the actually substantial atomic binding energy (Sec. II B). 3). A 
mixture of two fluids interpenetrating on an atomic scale lacks the mechanism for the main- 
tenance of the atomic bonding characters of the respective fluids (Sec. II B). 4). The critical 
velocity criteria for superfluidity of Landau^^ and Feynman"^^ (1955) have both incorrectly 
ignored the primary low energy excitations due to phonons (which in fact sensibly presents in 
the LLBF superfluid dynamic picture), and the comparison of velocities of different masses 
they made does not conserve energy (Ref. 55). 5). The "roton" picture does not pre- 
dict the experimentally observed overall features in the excitation at (Sec. II B). 6). 
The "Bose-Einstein condensation" picture is itself not capable of explaining superfluidity.^^ 
These difficulties inherent in the foundation(s) of the LLBF theory, cannot to us just be 
removed by a few minor improvements. Therefore, as the motivation for beginning this 
work, we felt it necessary to develop a theory that starts from an alternative - an essentially 
consistent - theoretical grounding. 

To provide an overall satisfactory interpretation of the complex behavior of He II, we 
flrst, in Sees. II-III, carry out a systematic derivation of a (microscopic) theory by using, 
to the broadest possible extent, the key experimental properties of liquid ^He as input 
information. Three key questions are to be addressed: 1) do the superfluid atoms primarily 
perform localized oscillations, or, diffusion, or do they stand still (see Sec. II A) ? 2) whether 
and how the atoms are interacting (see Sees. II B, HI)? 3) when the two fluids coexist in He II 
in the higher T end, do the superfluid atoms aggregate in large regions such that they possess 
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a characteristic bonding energy, or, do they intermingle with the normal fluid on an atomic 
scale (Sec. II C)? The answers to these questions have been achieved in this work after 
making iterative contrasts between what is predicted by an assumed "theoretical model" 
and what is observed by experiment, and between the prediction of one property and its 
implications for others. The key components of the theory are therefore already characterized 
in Sees. II-III. In Sees. IV - VIII we establish the corresponding formal theoretical and 
quantitative description of He II, including the composition structure of He II (Sec. IV), the 
set-up of equations of motion of the superfluid atoms (Sec. V) and solving (VB-VC), and 
the prediction from the solution of the physical properties, including excitation spectrum 
(VB,VC), the stability condition of the superfluid (Sec. VI), the quantization of circulation 
(Sec. VII), and the scattering function near the superfluid bond excitation (Sec. VIII). 
The superfluidity (its mechanism and the critical velocity) is treated in a separate paper, 
11.^^ These above contents constitute the microscopic part of the theory of superfluid *He. 
Sections IX-X present a corresponding statistical thermodynamics description of He II. Brief 
description of the theory has previously been presented in unpublished internal reports'"'^ 
except for Ref. 57 (2000) which mainly discusses the experimental probe of the QCE effect 
as fully presented in paper II. 

II. THE FORMULATION OF THE MICROSCOPIC THEORY BASED ON 
OVERALL EXPERIMENTAL OBSERVATIONS 

The static structure factor, S{q), of the (pure) superfluid, as has been probed both by 
neutrons"''*'"''^ and by X-ray,^° exhibits a heavily damped flrst peak at about qi, « 1.93 
A^^and almost vanishing higher rank peaks (FIG. lb). Immediately, this signatures that the 
superfluid retains only a short range ordering, with a mean atomic separation a = 27r/g(, = 

3.3 A. The detail of this ordering is not our concern; for a large part of our purpose here 
it suffices to conveniently, fully represent the atomic configuration with three orthogonal 
arrays of atomic chains in X, Y, Z directions, we thus have a simple ciibic pseudo lattice with 
a lattice constant a and a pseudo lattice reciprocal lattice vector b = 2TT/a = qt; see further 
Sec. VB. 

Experiments, including specific heat (Sees. I, II ), frequency dispersion (Sec. I), and second 
sound, have commonly shown that He II is an coexistent of the superfluid between about 
2.17 to 0.6 (^1) °K, reduces to predominately a (pure) superfluid below 0.6 °K, and at 
zero °K is himdred percent a superfluid. Below in Sees. IIA-IIC we infer from experiments 
the excitation and dynamic schemes in the two respective temperature regions, and the 
microscopic composition structure of the two-fluid coexistent and the related dynamics. 

A. Thermal excitation via single phonons and atomic dynamics scheme of the 

superfluid equilibrium 

(i) The dynamic structure factor The most revealing indication as to the underly- 
ing dynamic scheme in the superfluid equilibrium, which predominates below < 0.6 °K, 

is from the dynamic structure factor S(q,uj). The measurements from inelastic neutron 
scattering^^'^''"''^ have shown that the S{q, ui) of He II is characteristically different at T 
below and above T\, at a given momentum transfer q, with q raging from up to at least 

2.4 A"^.''^ The strongest distinction presents at an intermediate q (~ 1.1 A^^).^*"'^" Here, 
the S{q,u)) exhibits a superposition of a sharp and a broad peak on the w axis, w being the 
neutron energy transfer. The broad peak, which has a non-vanishing intensity at w = 0, 
at T above T\ identifies with the broad diffuse peak of the normal fluid He I and reduces 
smoothly in intensity as T is reduced to below T\. The sharp peak appears abruptly as T 
falls just below Tx; its intensity increases as T falls and then basically stabilizes from about 
1 °K downwards. 
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FIG. 1. (a). The experimental static structure factor S{q) 60 is plotted to assist discussions 
in the text. (b). Total thermal excitation energy e{q) versus the neutron momentum transfer q 
(equal to phonon wave number K) for the superfluid of He II. Circles indicate data points from 
inelastic neutron scattering measurement at 1.12 °K67. The solid and dashed curves show the 
theoretical dispersion (Sec. V); see further FIG. 7. The dotted curve illustrates a regular (i.e. with 
no superfluid bond excitation and with no nearest-neighbour repulsion) longitudinal single phonon 
energy dispersion, given by Eq. (26a) for all K values, obtained with a = 3.3 A (27r/o = 1.93 A"'^). 
(c). Lower graph: ratio of neutron scattering intensities, /i(Ti) and 72(72), at two temperatures Ti 
(< Tx) and T2 (> 2a) as indicated in the upper graph. 



It is natural to correlate this T dependence of the scattering intensity with the two-fluid 
fraction variation just mentioned. Woods and Svensson^^ based on quantitative evaluation, 
have shown that indeed the sharp and the broad peaks are attributable to the superfluid and 
the normal fluid respectively. Now, for the superfluid of concern, we can thus concentrate 
oursevels below only on the sharp peak. The compelling feature with the sharp peak is 
that, with the observed finite width being in principle only due to the finite instrumental 
resolution, it has the three characteristics of the phonon excitations in a crystalline 
solid: it is qualitatively narrower, it has a high frequency cut-off and, it peaks at a finite to 
and has a vanishing intensity at zero oj. Such a sharp peak in S{q,uj) which represents the 
probability that an energy transfer of u occurs inevitably implies that, the excitation must 
have occurred between two discrete energy levels. Experiments using pulse transmission 
techniques'^® have on the other hand shown that, at long wavelengths, the excitations are 
via sound waves. The combined information above hence characterizes the excitation to be 
via the creation and annihilation of phonons, each having an energy quantum, £{q). We 
refer to these as "single phonons" to separate from the "collective phonons" of Sec. II B. 
It should be remarked that, phonon excitation, with a well defined atomic dynamic scheme 
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in Dcbyc's theory of solids (1912),^^ is terminological also held in the LLBF theory for He 
II but, as will be illustrated in FIG. 2 below, it clearly is unable to be produced by the 
dynamics scheme of the LLBF theory. 

To infer from the above feature of excitation the underlying atomic dynamic scheme, wc 
shall first be guided by a general principle of quantum mechanics, which states that the 
state of a quantum particle will be discretized when under spatial confinement. In a bulk 
superfluid a condensed atomic assembly, if and only if assuming the atoms are localized - in 
both ends of its excitation, then each atom is wholly confined by its surrounding atoms, as is 
in a solid. And only then, the atom will then acquire fully discretized states; the transition 
between these states will result in the sharp peak of the S{q,uj), as is how such a sharp 
S{q, u)) peak in a crystalline solid (at low T) regularly results. (Although, the two present a 
few important distinctions which we discuss in Sec. V B.) Based on the above arguments and 
furthermore on the variety of other uniform indications to be described below, we propose: 

In the superfluid equilibrium, the thermal excitation is predominantly via the creation 
and annihilation of (longitudinal) single phonons, whcih predominate at temperatures 
below 0.6 (~ 1) °K. This excitation can only have occurred within the (pure) superfluid 
and resulted from localized (harmonic) atomic oscillations (of an amplitude 
u(R)), referred to as the LAO dynamics scheme. Each oscillator roughly corresponds 
to a single helium atom, separated by the interatomic spacing of He II, a ~ ,3.3 A, and 
oscillates in a potential well V (of a depth of ~ —7.2 °K/atom as will be detailed in 
Sec. II C). 

Based on the LAO dynamics scheme, the solution of equations of motion of the superfluid 
atoms in Sec. V actually reproduces the experimentally observed excitation feature, i.e. the 
excitation is between discretized energy levels e„, n = 0, 1, . . . (FIG. 2a), each differing by 
an energy quantum of a phonon, £{K). In effect, LAO is the only dynamic scheme able to 
give rise to a sharp diffraction peak for the superfluid. In contrast, the LLBF picture of 
phonon excitation (FIG. 2b), namely the excitation results from the knocking out of atoms 
from the Bose-Einstein condensate into the normal fluid, will incorrectly predict a broad 
peak of the dynamic structure factor. Evidently, this is because the energy of an atom in 
the normal fluid end, E{p) at a given momentum p, and consequently the excitation energy 
E{p) — Ebec, are continuous, despite the single valued BEG energy Ebec- 




FIG. 2. Excitation scheme for the superfluid equilibrium as proposed: (a) by the present theory. 

The excitation is via the creation of a single phonon, corresponding to a transition between the 
discrete energy levels en{q) and e„+i(q) of the oscillation mode q of the localized (relative to ^) 
superfluid atoms, (b) by the LLBF theory. The excitation of the helium atomic motion is produced 
by a transition from the single degenerate state of the " Bose-Einstein condensate" (the superfluid) 
of energy Ebec, to the continuous normal-fluid energy E{p) at a given momentum p. 

The inference of localization and the LAO from the sharp S{q,uj) peak is further sup- 
ported by the observation of a one-to-one link between a sharp S{q,u)) (or broad) peak and 
atomic localization (or diffusion) which being general and system- independent, as has been 
identified, for the first instance, by Skold and Larsson''^ in their comparison of the scattering 
features from the solid and the classical liquid phase of argon, and similarly from the two 
different phases of other systems. Where, the broad scattering intensity at lower energies 
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was assigned to be due to difFusion, the so called " quasielastic scattering" . In another ex- 
ample, a sharp S{q,u) peak centered at a finite co for a fixed (intermediate) q has also been 
observed in various metallic liquids, such as lead by Egelstaff'^^ and rubidium by Copley et 
a.\7^ Copley and Lovesey*^*^ correlated that feature to the stronger metallic bonding of Rb 
compared to that of a classical liquid. We can readily infer that the strong bonding functions 
to confine - to localize - the atomic motion. 

The LAO scheme also immediately explains as to why mainly the intermediate q scattering 
rather than the low q is specific with the superfluid; at a g < 0.4 A^^, for T below and above 
T\, the position, intensity and width of the S{q,u)) peak of liquid helium exhibit no distinct 
changes except that the rate of increase of the width is slowed down.™ Since K {= q) ^ 1.1 
corresponds to a wavelength A (= 2TT/q) ~ 10 A crossing over only a few atomic spacings 
(with a ~ 3.3 A). Hence a coherent oscillation, if to occur at such a short wavelength, would 
require all of the atoms move cooperatively and the atoms must therefore be localized. By 
contrast, at long wavelength crossing over many atoms, the oscillation occurring by means of 
density fluetiiation is insensitive to diffusion of some individual atoms. Indeed, in addition 
to the observation for liquid helium above, neutron scattering experiments have shown that 
the S{q,Lj) functions from normal liquids all present a certain sharp peak at g < 0.4 A"^.^*^ 

(ii). The excitation spectrum The LAO scheme is further pointed to by the excitation 
spectrum (FIG. lb), a;(g) due to the sharp peak of S{q, w)^^ in two respects. First, except for 
the peculiar non-zero value, As ~8.6 °K at qt ~1.93 A~^, w(g) has a sinusoidal dispersion 
which is the solution of the equation of motion of (relatively) localized helium atoms, Eq. 
26 in Sec. VB, for the longitudinal mode. Second, the solution informs the linear sector 
to be due to propagating elastic sound waves, here phonon waves, of a velocity 238.9 m/sec 
which is in near exact agreement with the first sound velocity ci = 239 m/sec from pulse 
transmission experiments.^^ Moreover, the experiment by Woods showed that the slope here 
is independent on T,^^ and disproves the postulation^^ that the slope might be instead due to 
the two-fluid density variation with T {ci{ps/ p)^/^). Furthermore, the dispersion "period", 
qb = 27r/a ~ 1.93A~^, giving a ~ 3.3 informs the oscillator size Op ~ 3.6 A,^^ which is 
just a helium atom. This can be similarly seen from the upper limit, q ~ 0.6 A~^, of the 
propagating phonon wave which development would require at least be a few oscillators 
within a wavelength. In contrast, the part of the "in- phase" oscillation of a large (Nm) 
number of superfluid atoms (to be discussed in the next section) , which alone is similar to 
the weekly interacting BEC wave or the zero sound wave''° in the LLBF terminology, has a 
maximum dispersion period Km oc 2TT/NMa « qb and an excitation energy negligibly low 
compared to the one of the single phonon here. 

(Hi). The T^ specific heat The LAO scheme is indicated by the Debye T^ behavior of 
the specific heat Cy of He II between ~ 0.6 °K (FIG. 3a). Insofar as a non-conducting 
condensed atomic system is in question, a Debye T^- Cy is broadly known to be necessarily 
given rise to by an atomic assembly where the atoms are localized and perform harmonic 
oscillations; the superfluid is not an exception as Sec. IX D explicitly treats based on statis- 
tical thermodynamics. By contrast, a classical fluid consisting of diffusive atoms generally 
shows a Petit-Dulong heat capacity. It follows that, the LLBF excitation scheme will in- 
evitably yield a specific heat strongly bearing the feature of the normal fiuid, which cannot 
be of a T^ behavior. Furthermore, the theoretical Debye Cy(T) obtained in Sec. IX D, 
when fitted to the experimental Cy data, produces a longitudinal phonon velocity of 241 
m/sec. This value agrees (within experimental error) with the ci value (239 m/sec) from 
pulse transmission^^ as well as phonon velocity (238.9 m/sec) from the excitation spectrum 
discussed in (ii). 
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FIG. 3. Specific heat of He II, Cv{T). Circles represent experimental values (Refs. ??); solid 
lines are statistic thermodynamic results of Sec. IX as based on the present microscopic theory 
of Sees. II-V. (a) for T below 0.6 °K, the theoretical Cv(T) oc r-\ Eq. (74)', is given by the 
derivative of the internal energy obtained from a canonical ensemble average over the states, with 
a Debye type of density of states, of the independent, indistinguishable single phonons, which 
are produced due to in real space the harmonic oscillations in the normal coordinates of the (N) 
independent, localized suporfluid atoms relative to ^. (b) from 0.6 to just below 2.17 °K, the 
theoretical A-Cv(T) (Eq. (74)') results principally from the conversion of the superfluid to the 
normal fluid, via the excitations of superfluid bonds. The lower and upper limits of To, 2.17 and 
3.0°K, yield the uncertainty boundaries of the theoretical Cv{T) shown by the solid lines 1 and 2; 
To = 2.6 °K gives the optimum Cv(T), solid line 3. 

(iv) . The inferred diffusion rate The LAO scheme is directly supported by the sohd-hke 
low diffusion rate in He II, being 10~^ cm^/sec as Eq. 26 in Sec. VC shows. The associated 
diffusion barrier implies a negative potential well in which the atom executes oscillation. 

(v) . The capability of predicting superfluidity The preceding observations inform that 
the superfluid atoms are not standing- still (in either the absolute or relative coordinates), and 
rather, the superfluid is fully thermally excited, with the excitation states, as Sec. IX shows, 
exhibiting a Planck distribution. The LAO scheme and the implied property above are 
furthermore supported by the fact that, as paper II shows, on the basis of the LAO scheme 
the QCE superfluidity mechanism can consistently, quantitatively predict superfluidity and 
critical velocity both in satisfactory agreements with experimental data. 

The one-to-one connection between an excitation spectrum as FIG. lb for He II and a 
fully excited He II is an inevitable truth, if one compares it with the analogous excitation 
phenomena of the vast body of crystalline solids, in terms of the phonon feature in general 
and the Kci behaviour in the linear sector in particular. Had Landau's criterion of critical 
velocity for He II, namely that if Vs < ci there is no phonon excitation, reflected the truth, 
then one will be led to the inevitably untrue conclusion, that any piece of crystalline solid 
moving, say against the surface of the ground, with a speed less than its sound velocity 
would be in superfluidity motion. 
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B. The superfluid bond excitation and atomic dynamics scheme; the gradual phcise 

transition of the superfluid 

The atomic localization in the LAO dynamics scheme derived in Sec II immediately sug- 
gests that, in the superfluid the atomic bonding has undergone a qualitative enhancement 
compared to the van der Waals (vdW) " bonding" energy of the normal fluid He I. We repre- 
sent the former by superfluid bond, Ab, which is defined to be the negative of the binding 
energy per superfluid atom that is given rise to by the attraction from many superfluid atoms 
within their correlation regime, T. Clearly, is directly associated with the stability of 
the superfluid phase. We propose: 

A quasi equilibrium process is predominated by excitations of superfluid bond (A;,) in 
two closely related Ab dynamics schemes, where Ab represents an attractive energy 
barrier which prevents the reference atom 1) from difi^using into its neighbouring site 
within the superfluid when locally perturbed, due to thermal fluctuation or an external 
incidence such as a neutron, and, 2) from transforming into the normal fluid state when 
globally perturbed, e.g. by bulk heating-up between 0.6 ~2.17 °K. 

By the definition above, Af, can be expressed by 

Ab = -Uso- < e„ > +Us. (1) 

Where, Ugo is the ground-state cohesive energy (per atom) of the superfiuid, representing 
the depth of the potential well in which a He atom is confined (FIG. 4al); Ugo = —7.2 

°K/atom as given by thermal measurements at SVP.^*'"'^* < e„ > is the pcr-atom kinetic 
energy of superfluid atoms and is << Ugo as will be verified in Sec. IX. ««(< 0) represents 
the "long range" attraction energy by atoms in the entire fluid bulk, so an atom in (Wg, 0) 
can diffuse virtually freely. The superfluid atom is, as Sec. VA will show, localized with 
respect to a kinetic energy even up to —Uso, thus one can expect Us << Af,; hence we have 

Ab « -f/,0 = 7.2 (°K). (1)' 



(al) 

superfluid 
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FIG. 4. (al). The potential energy V{u) for a superfluid atom at Xi and for one separated at a 
distance a at Xi+i, solid parabolic curves, theoretically constructed as in Eq. (20)" illustrated for 
the mode Ka — 0.5n; see FIG. 8 for V{u) at different K values. In the potential wells, the superfluid 
atoms perform independent harmonic oscillations in the relative coordinates ^. The potential well 
depth Vo is given by the experimental ground state energy Uoa = —7.2 °K/atom 88. The activation 
of an atom from the bottom of the well to just above zero, at a energy cost As = At + Sy, represents 
the excitation of a superfluid bond A^, with a site creation energy 6v < e > represents the average 
kinetic energy per superfluid atom, and by Sec. IX < e ><< Uos, hence A^ ~ Uos- The dashed 
lines indicate the repulsive part of the potential, V{u) > 0, and is associated with the phonon 
excitation within 0.6 < K < 1.4 (1/A) where £{K) > At; cf. FIG. 1. (a2) The ground state 
wavefunctions (solid curves) resulting from the solution of equation of motion of superfluid atoms. 
The dotted curves show wavefunctions for Ar atomic mass. (b). Illustration of the normal fluid 
potential energy for two adjacent atoms. E = Un — Uon — 3.1 °K for T just below T\ is the kinetic 
energy, and Avdw the average Van der Waals bonding energy. The diffusion barrier height between 
the adjacent atoms, being about 1.6 °K relative to E, is estimated such that it produces a diffusion 
constant in the order of 10~* cm^/sec at 1 °K, a typical value for liquids. V = —Sen is the energy 
region in which the atoms are free atoms throughout the bulk. Uon = Uos is assumed. 

Cases (i)-(iii) below show further experimental indications of the superfluid bond. 
(i). The A specific heat The superfl.uid bond excitation upon a global perturbation is 
directly indicated by the experimentally observed A- Cy of He II between about 0.6 to 2.17 
°K (circles, FIG. 3b). For if the superfluid to normal fluid conversion, which occurs in 
large quantity here, energetically involves the conversion of to the qualitatively smaller 
vdW bond of the normal fluid (FIG. 4a,b), then inevitably a large excess of heat supply 
is required. The activation energy, denoted by A.^^, is thus consumed mainly for and, 
additionally for a site creation energy 6vn and the kinetic energy E of a, normal fluid He 
atom. Hence, 

A^„ = Afe + Syn + E. (2a) 

Using (2a), the thermodynamic evaluation in Sec. IX will justify the excitation scheme 
above to actually reproduce the experimentally observed X-Cy- oc exp(— ^^). 

Furthermore, the A^ excitation scheme is supported by the observation tSat, a qualitative 
bonding change when transforming from one phase to the other is the common underlying 
excitation scheme for the second order phase transitions of a variety of systems (e.g. alloy) 
with a A specific heat. 

(ii). The dynamic structure factor about qb The superfluid bond excitation upon a local 
perturbation is directly pointed to, in the first instance, by the peculiar features in the 
S{qb,uj), at qt ~ 1.93 k-\ of the superfluid (FIG. lb). 15.67,70,89,90 j^^^ of these are: 1) it 
peaks at a finite u), Lj{qb) = A^ ~ 8.6 °K which is basically independent of T. 2) its peak 
width is qualitatively narrow and on the T axis exhibits a sharp edge at Tx.^"^ By contrast, 
the peak of S{q, lo) of He I is broad, and its position diminishes towards zero as T rises, e.g. 
to ~4.5 °K/atom at T = 2.9 "K.^" On the other hand, a simple cubic crystalline solid, for 
example, would show a zero frequency at qt- The underlying Aft dynamics for the above is 
explained below. At qi,, which reciprocal being ~ 3.3 = a A as noted earlier, the solution for 
the equation of motion (Eq. 12) relative to ^, the dotted line in FIG. lb, gives no phonon 
excitation. And, in real space, as viewed by a neutron committing a momentum exchange 
of qb = Kb, all atoms are located at the bottoms of their potential wells; cf. Sec. VB(ii). 
Hence, when a He atom here is knocked by a neutron transferring an energy As just a 
fraction larger than the well depth \Uso\ ^ A(, (FIG. 4a), it is, if moved at all, knocked to 
just a fraction above the well. Then, the Ag is mainly consumed to excite the A^ and in 
addition, for a small site creation energy 5s- As similar to Eq. (2a), we therefore have: 

As = Ab + S,. (26) 

As the superfluid atoms are localized, —At, and Sg arc therefore relatively sharply defined 
except for a finite smearing, primarily in Sg , due to the short-range-ordering resultant fluc- 
tuations in the instantaneous atomic configuration. This therefore necessarily explains the 
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qualitatively narrower peak in S{q,uj). On the other hand, a finite smearing imphcd above 
is indeed supported by neutron scattering experiments which have shown the S{q, w) peak - 
though narrow - has a finite widtli (e.g. systematically presented in FIG. 8 of Ref. 67), being 
a couple of times the instrumental resolution as used e.g. in Rcf. 73, and being also at least 
several times wider than Landau's theoretical prediction based on roton-roton collisions. In 
contrast, in a normal fluid the essentially free atomic diffusion results in a significantly broad 
peak in the S{q,uj). 

The A;, dynamic scheme is further supported by the third feature in the S{qb,Lu), namely 
that its intensity is the highest compared to at other q values. Prom the experimental data of 
Refs. 70, 68, the relative intensity, i.e. the ratio of the intensities, Ii and I2, below and above 
Tx respectively is estimated to be 14 at qt (FIG. Ic), which is about two times the second 
highest intensity, /1//2 ~ 7, at g = 1.1 A~^. An amplification of the scattering amplitude at 
qb is to be expected, if an incident beam of monochromatic neutrons are indeed scattered by 
single atoms from an array of plans at a mean separation a, at which the condition = 27r/a 
implies an effective Bragg scattering (sec further Sec. VIII). 

That the same dynamics underlies both the A specific heat and the excitation at qb, as 
implied in (i) and (ii), can be justified by the fact that the latter indeed constitutes the 
major contribution to the former, as is shown both by our thermodynamic result. Sec. IX, 
and by the accurate numerical calculations of Yarnell^^ and others. It is noteworthy that, 
a frequency dispersion with such a minimum at qb was initially devised by Landau in his well 
known 1941 paper, in order to yield a A specific heat. We here see that, the correspondence 
between the Cy and the uj{q) does not necessarily imply "which" dynamics. 

(Hi). The ion mobility The Ab dynamics scheme upon local perturbation is indicated, 
in the second instance, by the experimentally observed behaviour of the He ion mobility n 
in He 11,93-95 .(^^r^ich obeys 

/i = Vi/\E\a (X exp(A/fcBT), (3) 

Ea being the applied field. Since by definition mobility is inversely proportional to the ion 
scattering rate and hence to the excitation probability P oc exp(— As/fc^T) (Sec. VD): 
namely it is oc 1/P oc exp(As/fcBr), which is just as represented by Eq. (3). Hence the 
A of (3) identifies with the excitation energy, which is determined in the experiments to be 
8.1 °K (for negative ions) or 8.8 °K (for positive ions). Vi is the ion drift velocity, having 
a magnitude typically 10 ^ 100 m/sec. Clearly, the ion has a much higher kinetic energy 
than the thermal drift motion of a helium atom which velocity is in the order of 10"^*^ 
m/sec (see Sec. VC). Hence in the ion drifting process, the ion is essentially not subject 
to a many-quantum-atom correlation; and also the ion principally will transfer its energy 
and momentum to the superfluid, similar to the thermal neutron scattering. Hence, as 
discussed in (ii) , the scattering is due to the Ab excitation. Indeed, the basically equality of 
A 8.1 ~ 8.6 °K) from the ion scattering and the Ag of neutron scattering (Sec IIB(ii)) 
further supports the A^ dynamics scheme. 

The excitations in the two environments discussed in (i) and (ii)-(iii) can be generalized 
to be 

A,„ ~ A6 + 5„+£;„, (2) 

a omitted for the local perturbation, a = s, and n for the local and the massive bond 
break-up, Ea =< e >— is the thermal kinetic energy of a He atom in the superfluid, and 
Ea= Em the normal fluid of He II. Using in (2b) the known values of Ag and A;, and with 
~ gives 5s = 8.6 - 7.2 = 1.4 °K. 



C. The composition structure of He II, the momenton wave and the collective 

phonons 

The realization of the relative localization of the superfluid atoms (Sec. II A), and the 
formation and preservation of their large bonding (to be justified in Sec. II B) and its 



11 



many-atom correlation origin (Sec. Ill), all point to a composition structure of the two-fluid 
coexistent as described with a momenton model (FIG. 5): 

A momenton is firstly an aggregation of a large number of superfluid atoms, of an 
irregular shape and size of macroscopic scale. The motions of the atoms in it are pri- 
marily the u(R)- oscillations as described in Sec. II A, which are about sites R's fixed 
in the center-of-mass coordinates of the momenton, ^, here. The two-fluid coexistent 
of He II consists of tremendously many momentous adjacent to the normal fluid re- 
gions, or within large (pure) superfluid regions. Dynamically as a whole, in a He II in 
thermal equilibrium at rest, each momenton oscillates slowly with an amplitude ^(H), 
about its equilibrium position, H; this motion is equivalent to the simultaneous and 
in-phase oscillations of all of the atoms in it, which (being of a small proportion) are 
superposed to the prominent coherent, not-in-phase u oscillations of the single atoms 
of Sec. II A. Under a pressure difference, S describes the translation of a momenton, 
and hence the steady flow motion in the absence of convection. Finally, the excitation 
of a momenton occurs via the creation/annihilation of collective phonons, giving rise 
to the momenton waves. 



(a) 

momenton (a superfluid aggregate) 




(b) 




FIG. 5. (a). A snapshot of the schematic composition structure of He II at a T at which both 
fluids have a significant fraction. Each superfiuid aggregate, a dotted region, represents a momenton; 
the momentons are randomly distributed through a normal fluid background, (b) . In a momenton 
region, out of the thick framed region in (a), the absolute coordinates r of a superfluid atom is 
represented by the vector sum of four relative coordinates: u describes the atomic oscillation relative 
to R, R the translation of the equilibrium position of the atom viewed from ^, the coordinates of 
the center-of-mass of the momenton, ^ the oscillation of the momenton relative to H, and H the 
translation of the equilibrium position of the momenton. 
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In contrast, a mixture of two fluids which are intermingled on an atomic scale, as is 
assumed in all of the existing models (or, the LLBF theory), lacks the mechanism for the 
preservation of the interaction of the respective fluids. The formation of momentons also is 
energetically favourable in two ways: 1) given the largo supcrfluid atomic binding originates 
from a many-quantum-atom correlation (Sec. Ill), the gathering together of the superfluid 
atoms will allow an optimum atomic attraction, thus yielding the lowest energy of the 
system. 2) it requires less energy to turn a normal fluid atom into a supcrfluid atom, if 
the former is at the edge of the existing superfluid region rather than inside the normal 
fluid. Furthermore, based on the momenton model of composite structure, the two-fluid 
fractions can be quantitatively predicted with only the experimental density of liquid helium 
as an input data (Sec. IV), in satisfactory agreement with the statistical thermodynamics 
treatment in Sec. IX. 

Direct experimental indications of the momenton wave include: (i) A slow propagation 
of the elastic wave - the so-called second sound wave - has been observed, e.g. by Peshkov 
et al.^'"' and Atkins ct al.,''^ in thermal pulse transmission measurements of He II. It occurs 
only at low frequency (less than 1000 cycle/sec); it has a slow propagation velocity, C2 — 
20 m/sec within 2.17 ~ 1 °K, and from about 1 °K downwards, C2 increases rapidly and 
the pulse becomes substantially broadened. Despite the broadening, a finite width of the 
pulse is observed even with T being as low as 0.7 °K, which tends to suggest a finite size 
of the momenton in the superfluid. (ii). The thermal conductivity in He 11^^ is observed 
to be abnormally large and not defined by a temperature gradient. The momenton wave 
essentially describes the phenomena of (i)-(ii). For instance, the simultaneous and in-phase 
oscillation of many atoms (of a momenton) in eflfect will produce the instant propagation of 
a heat pulse; in other words, the pulse is at once carried forward by a distance having the 
dimension of a momenton. 



III. THE ORIGIN OF THE SUPERFLUID BOND 

The experimental evidence of Sees. IIA-IIB together points to the cooccurrence of the 
large atomic bonding and the (relative) localization of the atoms in the superfluid. Such a 
cooccurrence commonly takes place in such processes as liquid solidification, for which the 
physics is generally understood. On the other hand, the physics for the superfluid atomic 
localization, which occurs uinquely as a result of the transformation from the classical state 
He I to the quantum fluid He II, appears specific. From the basic fact that the superfluid 
atoms are quantum mechanical, we can infer their enhanced bonding mechanism as follows. 
A quantum (fluid) atom is marked by a broadly spread wavefunction tp{r), through a volume 
at least a few times greater than a (see Sec. V). In this volume, at any location we will 
find a fractional density of atom A, ip\{rAi) and that of atom B, tp'^^{rBj)] for simplicity, we 
call them partial-atom Ai and partial-atom Bj. Each partial atom consists of a fractional 
nucleus (of a charge ilP{r)Z'^^) and a fractional electron cloud (of a charge ilP{r)Z'^~) in just 
the way as the complete He atom, except for each being weighted by the fractional density 
'4P{r) {Jy tjj'^{r)dr = 1.) The above yield three characteristics in the atomic interaction. 1). 
the exchange-correlation is amongst all of the partial-electrons of all of the partial-atoms, 
as well as among all of the partial nuclei, that tend to approach the same locations. In this 
complicated correlation process, the surrounding atoms extending to many neighbors are 
strongly involved with one another. 2). given an AB atom pair (meaning A sitting to the 
left of B), some of the partial-atom pairs can assume the BjAi configurations. These BjAi 
pairs will repel with each other when A and B are tending to part. Both of the above aspects 
will lead to the observation that the atoms are prohibited from separating from one another 
and arc thus driven into a (relative) localization. 3). the interaction strength between a pair 
of partial-atoms is only a fraction of that of the total AB interaction. This will particularly 
affect, actually soften, the short-range repulsion. 

In addition, the localization allows a cooperative motion of many helium atoms, the LAO 
scheme. The atomic polarizations are therefore induced systematically, which essentially 
corresponds to the London or dispersion energy,^'''^* of the form e^huo/2a'^r^ , rather than 
randomly as with the normal fluid van der Waals interaction. 
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IV. QUANTITATIVE DESCRIPTION OF THE COMPOSITION STRUCTURE 
OF He II, THE TWO FLUID FRACTIONS 



Consider a He II bulk at a given T of the composite structure as depicted in Sec. II C; let 
the ith superfluid (s) or normal fluid (n) region in it be denoted by Vsi or Vni- Then one 
expects that the s and n component each has its own specific density, denoted by ps and 
Pn- Incorporating further the factor that the atoms of the normal fluid arc classical and of 
the superfluid are quantum mechanical, the local density of He II, i.e. the number of atoms 
per unit volume at r is thus 

or, 

p(r, T) = Ps(y, T)5{v e Vs^) + p„(r, T)5{v e Vm). (4)' 

where tp{rj) and S{rj) are the wavefunctions of a quantum and a normal helium atom re- 
spectively. The J2j '^uiis ^^^^ atoms within the wavefunction overlap region including 
= I". 

We can define the volume fractions of the s and n fluids respectively by 

/»(T) = 1-/,(T) = ^. (6) 

where the sums cover all of the regions of the respective fluids in the container, and V , Vs, 
and Vn are the total volumes of the two-fluid mixture, the s and the n fluid, respectively. 
The average density of the two-fluid mixture is therefore 

p(T)= (l/V) [ p(r,T)dr 
Jv 

= f,{T)p,{T) + UT)p^{T), (7) 
where, the average densities of the respective fluids are 



PiiT) = ^ = {l/Vt) ( pt 
Vt Jvt 



{r,T)dr, t = s,n. (8) 



Based on Eqs. (5)-(7) we now determine fs{T) and fn{T) using only the experimental 
densities of liquid ^He®^ as input data. Based on the close resemblance of the superfluid 

to a solid in a few important ways, e.g. as is shown by the essentially T independent 
thermodynamic functions (Sec. IX), we can assume that ps is basically independent of T 
on the scale of concern, and can then write ps ~ p{0), p{0) (= 0.1450 g/cm ) being the He 
II density at zero Kelvin. To obtain the pn of He II, wc next extrapolate Pn{T) of the He I 
phase into the He II phase region using the thermal expansion function: 

.^x 0.1495 

PniT) = o. (9) 

(1 -I- 0.001603T - 0.002035T2 0.001219T3)'^ 

Where, the numerical coefiicients are obtained from a least squares fit (solid line, FIG. 6a) 
to the experimental p{T) data of He I (circles) between 2.6 to 4.4 °K. Substituting into Eq. 
(7) with the experimental p{T) and the extrapolated Pn{T) given in Eq. (9) for T < 2.17 
°K, and by further combining with Eqs. (5)-(6), we then obtain for He II 

MT) = {p{T) - pn{T)) I (p, - p„(T)) (10a) 

f.{rr\ 0.1495 

_ ' (H-0.001603T-0.002035T2+0.001219T3)^ 

~ Qi^g 0.1495 

(1-|-0.001603T-0.002035T2+0.001219T3)^ 

and fn{T) = l- fs{T), (106) 
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FIG. 6. (a). The density pn{T) of the normal fluid component of He 11, soUd line within T <T\, 
as is extrapolated using the thermal expansion function Eq. (9). The coefficients are determined 
by a least squares fit to the experimental data of pniT) (circles) 99 between 2.6 - 4.4 °K. (b). The 
volume fractions /s (T) and /„ (T) of the superfluid and the normal fluid of He H, dashed curves, 
as given by Eq. (10), Sec. IV, established based on the composite structure of He II, namely largo 
separate aggregates of the two fluids, and with the experimental density of He II £is input. The 
number fractions Xs and Xn of the superfluid and the normal fluid of He II, solid lines, obtained 
from thermodynamic evaluation using Eqs. (70)-(70)', Sec. IX. The lower and upper limits, solid 
lines 1 and 2, are obtained using To = 2.17 and 3.0 °K, and the optimum values, solid lines in the 
upper graph and solid lines-3 in the lower graph of (b), are obtained using To = 2.6 °K. 

where all densities are in g/cw?. The resulting /^(r) and /„(r) (dashed lines in FIG. 6b) 
vary rapidly with T from about 50 % to 100 % and conversely between 2.17 - 0.6 °K, and 
approach to constant below 0.6 °K. In comparison, by Eq. (9) Pn{T) and ps{T) vary by 
less than 4 % across the entire T range of He II. Combining the above features with (7), 
we therefore obtain that the T-dependence of p{T) throughout the transition in 0.6-2.17 
°K is principally determined by the variations in fs{T) and fn{T) with T. It is further 
particularly noticeable that at T\, fg (^ 0.5) and /„ (^ 0.5) are not zero and one. The 
approximations involved would not cause any significant error as has been tested; rather, 
the non-unity of the fractions at the Tx is plausibly as a consequence of the second order 
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phase transition. For, in such a transition a measurable n to s conversion begins to occur 
at a T, denoted by To, well above Ty, that is, fs{To) = and fniTo) = 1. The physics 
involved will be further clarified in Sec. X. The fs{T) and /„(T) obtained above are in good 
agreement with the thermodynamic results for Xs,Xn in Sec. IX, tlic! resulting specific heat 
and entropy obtained based on which are in good agreement with experimental data; and 
these are consistent with the characteristic variation across the He II temperatures exhibited 
by a variety of properties of He II, including specific heat, and second sound. 

V. EQUATIONS OF MOTION, THE SOLUTIONS, AND THEIR USE FOR 
DIRECTLY PREDICTING PROPERTIES OF HE II 



The experimental evidence discussed in the preceding sections has pointed to the separate 
manifestations of the superfluid atomic dynamics with differing time and energy scales. 
Namely, a momenton oscillation, being of the scale of ~ 10^ cycles/sec, is negligibly slower 
compared to the single atomic oscillation with respect to ^, typically of the scale of 10^^ 
cycles/sec, so that the former can be regarded as being stationary within the oscillation 
cycle of the latter. Furthermore, the diffusion rate of the helium II atoms, with a diffusion 
constant ~ 10~* cm^/sec. (cf. Sec. VC), is negligibly slower compared to the single atom 
oscillation rate, and will on the other hand neither effect the momenton oscillation. We now 
establish the equations of motions accordingly via decomposing the total motion in a similar 
way to the above. First, we may express the absolute position r of an atom with the vector 
sum 



where u is the displacement of a He atom relative to R, R is the equilibrium position of the 
atom with respect to ^, ^ is the displacement of the center-of-mass of a momenton relative 
to H, and H is the equilibrium position of the momenton with respect to the laboratory 
system (FIG. 5b). Accordingly, the total interatomic potential is V"(R + u) + V(H + ^), 
with the 1st and 2nd terms being due to the atoms, as viewed from ^, and the momentons. 
As discussed above, the relative coordinates are effectively decoupled, meaning dr/dx = 
du/dx + dR/dx + d^/dx + da/dx, etc. Hence, we obtain four separate equations of motions 
for an atom at site I of mass m and an momenton j of mass Mj, which in the semi-classical 
form formally are the Newton equations: 



A. The equations of motion via coordinate decomposition 




(11) 



Fr = m 



Fi=m 



d^ni{K,t) 



(13a) 



(12a) 



/ = 1,2,... 




d^i^{KM,t) 



(14a) 



(15a) 



j = l,2,.... 



Quantum mechanically the total Hamiltonian and the wavefunction are 



H = Hu + Hr + H^ + He, 
^(r) =^(u)</.(R)fj($)!Z'E(H). 



(16) 
(17) 
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Here Hu, Hn, and -ffs, and (/?,(/), VPj and are the Hamiltonians and wavcfunctions in 

the respective coordinates, u, R, 4 and H. The counterparts of Eqs. (12) - (15) are the 
Schrodinger equations: 

Hu^ = e{K)^ (126) 

HR<i) =c {K)(t> (13&) 

H^$^ =e {KM)^i (146) 

FH^H(S) = e,(fc,)#E(S). (156) 

In Sees. V B-V E we set up explicitly for the above equations and solve, and based on the 
solutions predict a few immediate properties of the system. 



B. Cooperative oscillation of single atoms and single phonon excitation 

1. Set-up of Eqs. (12a)- (12b) 

Consider in the ^ coordinates the supcrfluid atoms are, as discussed in Sec. II A, localized 
and each oscillates about its equilibrium site R(X, Y, Z) along a chain of atoms say in the 
X direction, in a potential well V{X) (FIG. 4a). The mean separation of the atoms is a; 
as based on the m, (Sec. II), a = 27r/gf, =3.3 A. The complete oscillation motion can be 
projected along three orthogonal chains in X, Y, Z directions, which can be oriented in an 
arbitrary direction as the fluid is isotropic. In three dimension the dynamics above can be 
accordingly represented for a simple cubic pseudo lattice (see also Sec. II), which we shall 
use below, of an apparent and effective lattice constants, a and (to be explained below), 
respectively, with one atom per unit cell. When needed, the effect of fluctuation about 
a may be plugged in e.g. by damping the average result obtained this way. As the fluid 
has no shear elasticity, the different parallel atomic chains are not correlated in motion or 
coordination; and along each chain only longitudinal waves present. 

The exact many-body potential energy of an atom at Xi is V{Xi-\-ui), which is a function 
of the displacements of all (A'"x) atoms within a correlation length T: 

V{Xi+Ul{uQ) = V{Xi -\- Ul{ui.^);Ui^,U2s, . . . ,UNr,s) (18) 

where ui^ = J^i' [{-^i + ui) — (A;/ + u/')] is the displacement of the atom at Xi relative to 
all other atoms in the chain. V{Xi + w;) expands into the Taylor series: 



dV 1 d^V_ 



where Vq is just the semi-classical ground-state potential energy with a zero-point energy 
effectively included; Vq will not contribute to the force except when the motion of AT; is in 
question in which case it together with the A';(t) motion is then separately represented in 
Eq. (13a) In the LOD scheme, ^^w^j, will average to zero for the atom oscillating about its 

equilibrium site. Because the superfluid bonding (Sec. Ill) has a large spatial extension, we 
expect that at larger ui values, varies only slowly with -uj, hence higher order terms than 

d'^V/duf^ are small altogether and can be neglected even at a relatively large ui. V{Xi +ui) 
then reduces to the quadratic 

n«0 = K) + ^^<=l^o + ^ai<; (20) 

where A;, being inexpHcit, is dropped for simphcity, ui = ui{ui^) and ai = 4-^- For the 
many-body representation (18): 
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_ (fV_ _ f dV 

= E ^ I ^ ' (21) 
where = = with all other terms cancelled out as the displacements 

is is 

altogether are symmetric about Xi. Taking derivative with respect to ui on both sides of 
(20) we get the force acting on the atom at Xi: 

- - V— ^^^^^ ~ -a u; (22) 

Eq. (22) or (21) contains only the non-zero, diagonal terms q = q' {q, q' = X, Y, or Z) of 
the block matrix components — — , representing the absence of shear elasticity of the 

fluid. Further we sum for ui^ over only nearest-neighbors: ui^ = —ui+i — ui-i + 2ui. This 
together with the last expression of (22) represents an effective nearest-neighbor interaction 
and is used below; and the many-body correlation effect will be effectively included in an 
effective ai (or ci) value determined by fitting to experimental data. With the ui^ above in 
(22), and in turn (22) for Fi in Eq. 12a we get: 

-ai{-ui+i - ui-i +2ui) = m-^, 1 = 0,1,... (12a)' 

this being the independent harmonic oscillator description for the cooperative motion of the 
superfluid atoms in the ^ coordinates. It is noteworthy that, first, the averaging over many 
neighboring atoms in Eq. (21) results in an identical ai for all sites at all time, despite the 
reality that the physical parameters at each site in general fluctuate. This thereby leads to 
the well-defined, site-independent excitation energy as from Eq. (12a)' (similarly Eq. (12b)). 
Second. Fi is a restoring force, either attraction or repulsion, as long as sign(F;) = — sign 
(u/j,), and further for the given |F;| oc \ui^\ here, V{ui) = — J Fidui^ will then yield Eq. 
(20); we shall encounter this notion below. 

From the above and using in advance the wjj, expression from solving Eq. (12a) in Sec. 
VBc, we have V — Vq = ^uf = ^muj'^uf, where ai(K) = m,uj'^{K), and can then express 

2 

for the total energy: Uk = + {V — Vq) for the normal modes of each oscillator. The Uki 
with pi = —ih-^, directly leads to with which in Eq. (12b) we have: 

2. Solution of Eq. (12b), single phonons and their wavefunctions 
Equation (12b)', which we here discuss first, has the well-known eigen energy solution: 
en{K)= (n+l-]nu:{K) n = 0,1,2,..., (23) 



2^ 

which consists of n phonons, each having an energy quantum 

£{K) = nuj{K). (24) 
The stationary eigen functions are (here u = ui) 



Vn{u) = CnHn [ \ j -j^ u \ exp( — —u ). (25) 
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Hn are the Hermit polynomials and c„ the normalization coefficients. The expectation value 
< e„ >= / $;(u,t)e„$„(u,f)du = J^nY^K^niK), where $„(w,t) = ifi{u) exp{iH/nt) are 
the time dependent eigen functions, yields the total energy of the system, which will be 

given in Sec. IX by an ensemble average, the internal energy Us — Uso- 

Equation (25) informs: 1) The spatial extension of ip{u), as can be measured by A„ = 

2w|,^(m)=o.5<^(o) = 2^ ^ ln'2, is large for the small hehum mass. At a low energy, e.g. 
u) = 1 °K, A„^^ ~ 8.2 A and is a few times larger than a. The minimum-energy (of the 
u motion) wavefunction ipo{u) — ^/^^e"^" at this uj is seen to extend into and beyond 
the first nearest-neighbor (solid curves, FIG. 4a2); and similarly with the excited states. By 
contrast, the wavefunction for an argon atom, of a mass mAr ~ lOmne, (dotted curves in the 
figure) has A„^^ ~ 2.6 A, which is less than the interatomic spacing OAr — 5.3 A taken for 
solid argon. A„, being a similar measure to the thermal de Broglie wavelength, thus defines 
that, (viewed now in ^) a helium is quantum mechanical and an argon classical. A„ increases 
with cj, but by Sec. VB2 at a larger lu the atom is then confined in a narrower potential 
well; and at a high co comparable to A;, the atom converts to the state of a normal fluid 
atom. 2) It follows from the large A„ that the single particle momentum of helium cannot be 
determined better than ?i/A„. 3) The immediate consequence of the large A^, which causes 
the wavefunctions of several nearest-neighbor atoms or so to be entangled, is to result in the 
many-quantum-atom correlation scheme, and the consequent substantial superfluid bond 
and the localization of the He atoms, as discussed in Sec. III. 4) However, A„ is only several 
times a or so say within a region T; hence the helium atoms are "indistinguishable" effectively 
within only the limited T region (T < or << y of the container), in which the atomic 
wavefunctions overlap. This is in contrast to a free quantum particle whose uncertainty 
region consists in the entire volume of the particle assembly (the free electron gas is an 
example of such). Furthermore, since the atoms arc localized and each can most probably 
be found about its equilibrium site, thus the individual atoms are not "identical" in respect of 
the sites they can occupy. Restating this condition, by regarding the localized atoms as being 
effectively " identical" instead, wc have that the centres of mass of the atomic wave packets 
are effectively "distinguishable"; the atoms will become "effectively identical" if one views 
each " atom wavefunction -|- site" ; the " independent harmonic oscillator" representation in 
this section is in effect a formal transformation to this status. Compare the discussion 
relevant to the R motion in Sec. V C. 



3. Solution of Eq. (12a), the dispersion curve for q < 1.93 A ^, and the potential energy 

Equation (12a) can be readily solved to give the oscillation amplitude ui{K) = 
j^^-i{laaK+ojt) ^ thus ui-^ = 4u; sin^(-^^^), and the normal-mode angular frequency and the 
£{K) by combining with (26) (both in J) 

/ j^. o • ( 

UJ(K) = X — 2sm 

^ ^ y m \ 2 

and e{K) = £{K) = nw{K) 

where : 

a, K < 0.51, or 1.87 <K <Kb- f 1"^ 

(26a) 

Ur, 0.51 <K< 1.87 (266) 

For K < 0.5 A-i {Ka « 1), (26a) becomes 

£{K) ~ hciK. (26a)' 



The phase velocity here is £/{hK) = ci = yoi/m a, being equal to the propagation velocity 
of the sound wave, {l/h){d£/dK). Parameterizing Eq. (26a) using the pulse transmission 
value 239 m/sec^^ for ci, we get £{K) = 18.25Js: (°K) (solid line-1, FIG. 7a and also FIG. 
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lb) where K in A"'^, in satisfactory agreement with the nentron-probed excitation curve^'^ 
(circles) in this K region. The slope of the latter from a least squares fit is ci = 238.9 m/sec. 
The above result has been used as a fed-back information in Sees. II A-IIB. 




FIG. 7. (a) Theoretical total excitation energy e(K) versus phonon wavevector K of the super- 
fluid of He II, solid and dashed curves. The excitation at a < Kh — (j/2 is due to the creation 
of a single phonon, £{K) (solid curve 1), as predicted from the solution (26a-b) of the equation 
of motion of the He atoms performing independent harmonic oscillations relative to ^ in the po- 
tential wells V(u) shown in FIG. 4. The above K region further divides into K < 0.6 A^^(Eq. 

26a) and < 0.6/f < Kh — cr/2 A~^ (Eq. 26b), where the apparent parameter a = 3.3 A and 

the effective (ir = 2.64 A are, respectively, used as the interatomic spacing. In the vicinity of 
Kh — 1.93 A^"'^, i.e. \K — Kb\ < ct/2, £{K) (solid curve 2) primarily results from the excitation 
of a superfluid bond, of an activation energy As, superposed with a small single phonon excitation 
energy: £{K) = As{K) + £{K), where e{Kh) = As. Neutron diffraction data points at 1.12 °K®'' 
are represented by circles, (b) Potential enenrgy of a superfluid atom, V{u), as a function of K, 
Eq. (20)' with u — ar/2, plotted for several values from 2.5 to 3.3 A. At a fixed Ur, V{^ar) has 
a positive peak in a certain K region which shifts towards the higher K end with decreasing Ur- 
Only if ar « 2.64 A, the K region where V{'^ar) > (dashed curve) coincides with the K region 
where S{K) > As (dashed curve, FIG. 7a); this K region is bounded between (0.51, 1.87) A~^ as 
indicated by the two dotted vertical lines. 



For 0.5 < K < K}j{= qt) ^, ^{K) involves explicitly and here clq^ — CLj^ . O"!^ , as an 

effective quantity of a, was introduced when setting up Eq. (12a)' in order to now account 
for an apparent distortion of E{K) from the experimental curve if a = 3.3 A is used in Eq. 
(26). For we recognize that, a helium atomic oscillator, having a broad wavefunction as Sec. 
VBb showed, would overlap with its neighbours and hence repel each other at an effective 
border at I^a,,[ < \^a\ (FIG. 8). If the excitations indeed include oscillations of amplitudes 
larger than ia^, then the finite oscillator size, reflected by the repulsion length AA' in FIG. 
8, must be excluded from the apparent a. Thereby at these modes the atomic oscillator is 
effectively confined in (— ^a,., ia^)- In resulting in a £{K) (dashed line in FIG. 7a and also 
FIG. lb) in agreement with the experimental data in this K region, we find a^. « 0.8a = 2.64 
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A. We will shortly justify that with this value, V{u) indeed turns to be > at u = ^ar- 



V(u) (°K) V(u) (°K) 




FIG. 8. Potential energy V{u) versus m a of superfluid atom centered at Xi or Xi+i, Eq. (20)' 
witli Or = 2.64 A, plotted for several K values as indicated in the inset. For K = 0.51 A~^, from 
which towards higher K values £{K) > As (see FIG. 7), the potential turns to be repulsive, i.e. 
V{u) > 0, at the effective border u ± ^ (curve 3). a is the apparent interatomic distance of He II 
as given by the reciprocal of Kt: a = 2tt/Ki, = 3.3 A. 

Using in Eq. (20)' the ai = as determined from the linear K region relation and the 
ui^ expression, both just given, and with Vq — Uqs = —7.2 °K/atom as determined from 
experiment^* '^^ (note that by letting Vq — J/os the zero point energy is effectively included; 
see also Sees. II B, IX) we have 

= ^° + = -^-^ + (20)' 

where V{iLi) is in °K, ui in A, and K in A"-'^; o„r = 2.64 A for the parameterization; 
the resulting V(u) for selected K values are graphically shown in FIG. 8. For K < 0.51 or 
K > 1.87 A~\ V{u) < for all |u| values up to ^a^. For K = 0.51 (or symmetrically 1.87) 
or 0.51 < K < 1.87 A~-^, V{u) turns to be repulsive at u = ^a^ or earlier - we preassumed 
this in the preceding paragraph. Furthermore, the repulsive feature of V{u) with K in 
(0.51, 1.87) A~^ is also precisely as expected from the excitation scheme derived in Sec. II. 
Namely, since in this K region the excitation energy £{K) > Ag and £{K) corresponds to 
the max(|u|) for the given K, the £{K) here must therefore result from repulsion. 

Compute V{ui) as a function of K with ui = in Eq. (20)', for different values from 
2.5 to 3.3 A. We find that at a fixed a^, V{^ar) has a positive peak in a certain K region 
which shifts towards the higher K end with decreasing Ur (FIG. 7b). Only if ~ 2.64 
A, the K region where V{^ar) > (dashed curve) coincides with the K region where 
£{K) > As (dashed curve, FIG. 7a). Thus we have verified that the solution a^. 2.64 A is 
unique. 

Furthermore we note that at K = Kf, = 1.93 A~^, or Ka = 27r, all atoms move in phase: 
= exp^'f+i^2'7r"hLt) = ^^'^ = 0- That is, in the Kf, mode, in real space all atoms 
are at the bottoms of their potential wells as viewed, for example, by a neutron scattered 
due to this mode. The oscillations at the other normal modes would cause, with respect to 
the neutron, a snapshot of somewhat irregularly displaced "equilibrium sites" of the atoms 
which will add to their already irregularity due to a short range ordering. 

Finally, the character of the localized oscillation of the helium atoms, being clearly distinct 
from that in solids, has up to now been represented in the formal treatments and are re- 
emphases as follows: it is relative to ^ rather than H, it is in cages fluctuating in location 
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and dimension, and it is in a shallow potential well which switches to repulsion at larger 
oscillation amplitudes. 

C. Single atom translation and excitation of the superfluid bond 

1. The semi-classical representation Eq. (13a), excitation of the superfluid bond 

Equation (13a) describes, in a superfluid equilibrium in the absence of a pressure and 
composition gradient, the thermal fluctuation resultant translation (self-diffusion) of the 
equilibrium position R of a helium atom. This motion can be tracked by solving the simul- 
taneous equations of motion for N atoms with an interatomic potential Vq, e.g. numerically 
with the help of computer. Macroscopically, the motion of, say, N' tagged atoms form a 
thermal diffusion flux of a concentration c(R) ; the diffusion equation e.g. along X axis is 

D c(R, t) - = (the macroscopic form). (13a)' 

D, the diffusion constant, can be estimated by 

D= Ca^ — e-^''^^^ 
27r 

w 2.4 X 10"* cmVsec (26) 

where C is the coordination number, q-^sI^bT ^j^g Boltzmann factor and, by Sec. II B, 
As, =8.6 °K, is the activation energy of a superfluid bond. The numerical value of D is 
obtained using the previous values for a, and a;(= 10~® 1/sec), with ^ « 1 and T = 1 °K. 
This D value is as low as that of diffusion in a solid, informing that a superfluid atom is 
virtually localized. Accordingly the solutions of (13a) and (13a)', of which c(r) is known to 
be a Gaussian function, are not of interest here. D represents the probability for an helium 
atom to translate from its bonding state at R to that at a neighboring site at R' (which 
can statistically also be R), by overcoming an energy barrier height A(, + (FIG. 4a): 

P(R^ R') oc £) oc e"^. 

With He II, the P(R R') processes of practical interest are ones caused by external per- 
turbations, and in particular, the manifestations of the processes in the excitation spectrum 
which we discuss below, and in the specific heat which we discuss in Sec. IX. 

Upon a local perturbation the excitation ?A, K = Kf, is by Sec. IIB(ii) due to activation 
of the superfluid bond: 

s{K) = As =8.56 ("K) K = Kb (26c) 

And a,t a K apart from Kb but in a narrow vicinity < \K — Kf,\ < cr is, by combining with 
Sees. II A and VB, therefore 

e{K) = As{Kb) + £{K) « 8.56 + 18.18|iV' - 1.93| (°K), 

()<\K-Kb\<^ (26rf) 

Here, £{K) = ?ic\\K — Ki,\ is the single phonon excitation energy of Eq. (26) after the 
simplification: — sin(7r — ~ — (tt — ■^) for tt — << 1. For the parameterization: 
Ag = A;, -|-(5v = 8.6 °K and ci ~ 239 m/sec as given previously; K is in A~^. The resultant 
£{K) (solid line-2 in FIG. 7a and also FIG. lb) agrees fairly well with the experimental curve, 
except for its hard turn at as contrasted to the experimental smooth turn. This diffcrnecc 
is as expected, since the theoretical treatment does not include the effect of the fluctuation 
(resulting from a short range ordering) in the superfluid bond length which occurs in the real 
superfluid. Indeed, experimntal excitation data at high pressure (~25 atm)^°°, which causes 
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a "solidification" efi'ect, has shown a sharp turn at Ku similar to the theoretical description. 
(There are also shifts of the dispersion curve to a higher frequency and to a larger period, 
which can be understood similarly.) 

As K departs from (Kh — j , Kb + ^ ) to lower K values , the probability of the A5 excitation 
will drop rapidly, for the liquid density n(r) drops as r departs from R. So that e{K) = £{K), 
which is the purely phonon excitation energy as of Eqs. (26a)-(26b). 



2. The quantum mechanical representation Eq. (13b) 

The potential energy for the R motion for a reference atom at i? is: 

V„(<0) = {„^"-'^"- (27) 

Here Aq is the zero point kinetic energy of the R motion; Vq = —7.2 °K/atom as before. 
With Hn = — + Vo, Vo as expressed above, in (13b) we have in the X direction 

|l^+Vo)0=C^ |X|<| (136)' 



C6 \X\>^ 



The eigen function of an even parity is 



0(X) = Acos{KX), \X\ < I (28o) 
0'(X)=Ce-^'W, \X\>^ (28&) 



where hK = y^2rn(cr^jVo]y and K' = ^ . Applying the boundary condition at X = 
\a/2\, i^tan(^) = K'{a/2), which, when substituted with the expressions for K and K' above, 
writes as 77 = ctan(0 = , ^ , with ^ = K(a/2) and /? = -v/2m|Vo| h(a/2), the energy 

solution is given to be 



2m{al2frf 



C=- .^._' o 2 - (29) 



Starting with a guess value for Vq and thus for /3, put a = Ur = 2.64 A for the same reason 
as before, solve graphically for 77 and C, and thus C; use the resulting Co value in (31) 
(below) and Aq in turn in (27) to acquire a new Vq value; and execute this iteratively until 
Vo reaches constancy. We find, for Vq being small here there is only one, hence the ground 
state solution: 

Vo = -10.65 (°K), thus 13 = 1.749, 

770 = 0.68, Co = 0.3l7r. (30) 

With the 770 and Co values above in (29) we get the ground state eigen energy of the R 
motion: 

Co= -7.13 "K/atom (29)' 

With the Co we further get: 

Ao =Co -Vo = 3.45 °K/atom (31) 

The solution of odd parity is discarded as it does not fulfil the requirement of a symmetric 
(ground state) wave function for identical boson particles (Sec. VB-2); see further Sec. VF. 
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D. Momenton oscillation and excitation of collective phonons 



1. Set-up of Eqs. (14a)-(14b) and the solutions 

Equations (14) arc formally set up for every momenton masses, Mj, j = 1,2,..., with the 
inter- momenton potential V(H + ^), in the independent harmonic oscillator approximation. 
Since the momenton contribution to the total energy and excitation is negligible (see below), 
we will not solve its exact motion here and will only look at its long wave propagation, using 
an average mass M and with an effective force constant a2- Equations (14) thus reduce to 
a single equation similar to Eq. (12a)'; its energy solution is thus: 

£2{Km) « nc2KM, (32) 

where the wave number and velocity of a collective phonon are Km = 27r/Yl, Am being the 
momenton wavelength, and 



respectively. Apparently C2 corresponds to the velocity of second sound discussed in Sec. 
II C; a2, M = {NM fm, Lm = Nmu and Nm are the restoring force constant, mass, dimen- 
sion and number of helium atoms in one dimension of a momenton oscillator, respectively. 
a2i — a2/N^ represents then the effective force constant between two He atoms at the 
interface of two momentous. C2(T) is dependent of T via: 

dc2{T) _ dc2 ONm dc2 da2i , . 

dT ~ ONm dT da2i dT ' ^ ' 

Suppose at a higher T a momenton will most likely find itself immediately surrounded by 

normal fluid regions, the reduction of T is thus mainly to increase the number of momentons. 
Then Nm and a2i and, by Eq. 34 hence C2(T), are basically constant. Meanwhile, the 
interaction of a momenton with the surrounding normal fluid gives: a2i < cn- Immediately 
this gives: C2 < ci. As T is reduced (to about 1 °K) such that with a falling T, each 
existing momenton increases in both Lm and in its adjacency with other momentons for 
which a2i enhances, i.e. << and ^ << 0. Hence, by Eq. (34) ^ << 0. In 

other words, C2 increases rapidly with falling T, this being consistent with the experimental 
observation discussed Sec. II C. Furthermore, as with Eq. (23), the momenton excitation 
energies are £^ = (A/'+ 1/2)£m{Km), where A/" = 0, 1, . . ., and £m{Km) = hQ{KM) is the 
energy quantum of a collective phonon. Typically £m(Km) ~ 7.6 x 10~^ °K, as based on 
the upper-limit in the experimental second sound wave frequency, Q, 10^ cycle/sec. Hence, 
£m{Km) is negligible compared to the typical energy scale of a single phonon, hcu ^ 0.8 °K 
(~ 10^^ 1/sec). Accordingly the momenton oscillation amplitude |^{(0P would be basically 
constant on the scale of |'/?(w)P and can be dropped from (17). Finally, since Am is large, a 
momenton wave will be absent in a small sample. 



2. Momenton wave, the atomic dynamics of second sound wave and temperature fluctuation 

Consider a given normal-mode frequency ft of momenton oscillation that is equal to the 
normal-mode frequency w of single-atom oscillation; the latter consists of all wavelengths 
(much shorter than the momenton dimension) within the momenton. Two extreme cases can 
present. With some momentons, the n periodic process is in phase with the lu oscillation. 
Evidently, the ^l wave of amplitude An , which will superpose with the oj wave of amplitude 
A^j, will then increase the to A'^^ = A^^+Aq. And conversely with some other momentons. 
Therefore, in each of the former momenton regions the phonon population n and thus the 
temperature T are increased by small amounts n' — n (x 2^n^w and T' — T, as (^w)^ oc n, 
noting Aq « A^; and conversely in the latter regions. So that a temperature wave forms, 
which is one defining characteristic of the so-called second sound wave.^°^ 
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E. Superfluid flow, translation of momentons or of single atoms as of Eqs. (15)-(15b) 



Equation (15a) describes a steady flow motion of the superfluid in the absence of con- 
vection and is expressed for the translation of the equilibrium position H of a momenton. 
Under normal friction conditions, the flow is driven by an external pressure difference P1 — P2 
acting on a cross sectional area S of the flow, hence Fs — {Pi — P2)S. One portion of the 
is used to overcome the viscous resistance from the channel wall / acting on the flow 
through a flow- wall contact area Asy,; and the remaining to accelerate the flow inertia, the 
flow mass, which is 

M,=^Mi = N,m. (35) 

i 

Ng is the total number of helium atoms in the flow, and m - as before - the atomic mass of 
helium. Because < u >=< ^ >= and also < R >w constant (C) (<> denotes the time 
average of a variable over an infinitesimal flow translation), therefore 

<r>{t) = S{t) + C. (36) 

(36) implies that the translation of a single atom, or a momenton, or the flow, describes the 
same motion. If / « 0, which will occur when Vg is below a critical value (which criterion is 
given in paper II), Eq. (15a) then describes the non-dissipative superfluidity flow motion. 

is then entirely consumed to accelerate the flow velocity, for example at time t = t, to 
Vs{t) = {l/Mg) Jq [Jq{Pi — P2)S{t)dt\ dt. If from t > r, the driving pressure is removed: 
(Pi — P2)\t>T = 0- Then Eq. (15) and its solution become simply 

<S> <r> 

— dt2 — = — dt^ — = °' """"^ 

d <r> d <S> . . ^ ^ . 

Vs = 7: = r = Vs{t) = constant. (37) 

dt dt ^ ' ^ ' 

The classical equation (15a) validly describes the motion of the center of mass of an atomic 
wave packet. 

The translations of the quantum helium II atoms are more fully, or strictly, described by 
Eq. (15b). As ffs = p^/2m+ Vs and Vs = in the absence of non-uniform external field, 
the solution of (15b) is an "effective" planewave: 

^< r >) ^ Ce''-^^<''>-'^*\ (38) 

with "effective" referring to the fact that each single atom, if driven externally, is accelerated 
according to the relation (15), to a velocity Vg with an effective mass equal to the flow mass 
Mg. The wave number kg is connected with Vg by the de Broglie relation: 

kg = mvg/h; (39) 

kg is thus similarly effective. The translation kinetic energy of an atom is accordingly 

eg = mvl/2 = n^kl/2m. The flow energy is thus: Eg = ^Mgv"^ = NgCg = Ng^. It 
is noteworthy that, Eq. (39) holds only for the steady flow motion here, where periodical 
processes of each atomic waves are essentially linear, independent and hence obey the su- 
perposition principle of quantum mechanics. However, when some part (e.g. the interfacial 
layer) of the flow is being resisted as opposed the remaining part, the former will interact 
with the latter given that the superfluid atoms are strongly correlated. Then, each atom in 
respect of its flow motion is no longer a "free" particle. 

By knowing Vg and the density pg = l^l"^, with ^(r) as in Eq. (17) the superfluid current 
density is given accordingly 

is = - {vrm = PsVg. (40) 
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For the two-fluid coexistent of He II under thermal equilibrium, by Sec. II C each of the 
fluids are maintained in their own regions with no intermingling on an atomic scale, therefore 
each fluid maintains its specific density. For a steady, non-compressible flow, the continuity 
conditions follow accordingly 

V-j. + ^=0 and V-j„ + ^=0. (41) 

For a steady superfluid flow in a straight channel, the above becomes V • js = ^§f- = and 
V-j„ = = 0. For a steady flow of a two-fluid mixture maintained by a pressure difference 
there is Vs = Vn, and the island-like composite structure (FIG. 5a) should then maintain. If 
now the driving pressure is removed, so that, starting from the layers by the wall, = 
and only the (pure) superfluid will flow. A natural course via which this actually takes place 
is that the superfluid will break through the normal fluid and re-shape itself into narrow 
strips. 



F. The total and full excitation spectrum, the total wavefunction and total zero point 

energy 

a. The total and full excitation spectrum. The total excitation energy for a superfluid at 

rest is 

e{K)=S{K) + Ab{K) + 0{£M) 

~E{K) + Ab{K). (42) 

The full spectrum of £{K) in (0, 1.93) is combinatorially given by Eqs. 26a, b,c,d (solid 
curve 1, the dashed curve, and solid curve 2 in FIG. 7 and also in FIG. 1). The single 

phonon and the superfluid bond excitations arc, as the evaluation of specific heat in Sec. IX 
will show, characteristic in the T regions (0, 0.6) and (0.6, 2.17) °K respectively. 
h. The total zero point energy is: 

3o= eo + Ao + 0(eo) 

~ 7.5 + 3.45 = 10.95 °K/atom (43) 

where eo and Aq are given by (76) and (31). 

c. The total ground state potential energy is: 

$0= Vb - 3o = Vo - eo - O(eo) 

= -7.5 - 10.5 = -18.5 °K/atom (44) 

where Vo is given by Eq. (30). 

d. The total wavefunction. Of the V-'(r), Eq. (17), for an atom at r(.T,y, z), given the 
rapid u oscillation, say in the x direction, ip{u{X)) as of Eq. (25) will on the time scale 
of the R motion show as an time average to be an even function relative to X. \I'{(^) is 
even both for a similar reason and because l^'^P is essentially fiat (and shallow) within the 
region of wavefunction overlap. (t){X) of Eqs. (28a)- (28b) is even. Last, \1/2(H) is identical 
for all atoms in the (flow) bulk. Hence, overall, on the time scale of R motion, by which the 
atoms may exchange positions on activation, ■0(a;) is an even function with respect to the 
equilibrium position X\ and x can be an arbitrary direction relative to the isotropic fluid. 
Therefore, ^->{x) is spherically symmetric about R. The product of the single particle wave 
functions therefore yields a desirable, symmetric total wave function of N identical (in the 
sense of Sec. V B-2) particles of the fluid bulk: 

!Z/(ri,r2, . . . ,rAr) = ■!/'(?'i)^(?'2) • ■ • V'(''Af)- (45) 
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VI. THE STABILITY CONDITION OF THE SUPERFLUID 



From the supcrfluid bonding character we can infer: 1). The substantial superfluid bond, 
being larger than the van der Waals "bond" of He I (quantitative evaluation given in Sec. 
IX), implies a substantial internal viscosity, larger than that of He I. Thus the superfluid is 
not an assembly of non-interacting atoms and hence cannot be of a zero internal viscosity 
- as portrayed in London's BEC and in the unified LLBF theory as well. Meanwhile, the 
"weakly" (or sometimes referred to as "substantially") interacting superfluid atoms" and the 
"zero viscosity" which are simultaneously held in the LLBF theory, are two inevitable self- 
contradictory aspects. 2). Now that the viscosity of liquid He I can sustain a non-turbulent 
macroscopic rotation (below certain velocity) of He I which is a commonplace experimental 
observation, the larger internal viscosity of the superfluid He H can therefore certainly 
ensure a (more robust) non-turbulent macroscopic rotation of He II. The superfluid is thus 
not irrotational. 3). It therefore also follows that for the superfluid a zero total circulation, 
Ktot = 0, is not a necessary condition for non-turbulence; instead Ktot is finite and quantized, 
as has been experimentally demonstrated and to which a consistent theoretical grounding is 
provided in Sec. VII. 4). Accordingly, no local vorticity needs be imposed on the superfluid, 
as the Stoke's theorem, applying for a viscous-less fluid, does not apply to the superfluid. 
5). It also follows from points l)-3) that the grounds for proposing the "absence" of internal 
viscosity as an explanation for the non-dissipative superfluidity motion are lacking. (This is 
in addition to the fact that an internal viscosity does not necessarily imply which "flow- wall" 
viscosity, which we are however not to expand on here.) 



VII. THE QUANTIZATION OF CIRCULATION 

Consider a supcrfluid bulk in steady rotational flow motion with a tangential velocity . 
The trajectory of a flow atom < r >, Eq. (38), can be now conveniently projected onto 
the circular coordinate £ along the flow streamline, <> being here for simplicity omitted. 
The atomic wave, assuming being localized relative to the rotational frame and hence not to 
swapped off the stream line, in the circular loop, of a circumference L, must satisfy according 
to quantum mechanics the continuity condition: 

ksL = n2n. (46a) 

This then yields a circulation quantization as to be expressed in Eq. (49) below. The case 
here is seen to be equivalent to an electron wave orbiting about the nucleus, except for their 
different centripetal forces, 

Alternative to the above, to examine its microscopic detail we below look at the process 
from the standpoint of wave interference. We first determine an important parameter, the 
length of the rotational wave train of a helium atom, and this in theory is simply 

L^ = Nx- Xs (47) 

where is the wavelength and Nx is the number of As 's of the rotational wave train. The 
de Broglie relation (39) informs about Ag = 2'K/ks, and not the Nx. As to the physical 
acquisition of L^, let us think that the action of passing a Uks to a helium atom consists 
in passing it a field - in the form of a (purely) rotational wave train whose momentum 
equals hks ■ [We are tempted to think of this field being an electromagnetic wave, although 
this link would require a justification which goes beyond the context of this paper.] Since 
this wave is given as the solution of the Schrodinger equation (15b), we shall merely take 
this consequence. We can then think of the total atomic wave train now rotating along £ 
as the stationary thermal do Broglie wave train of a helium atom as being superposed in 
the rotational frame onto the purely rotational wave train of length Ag • Nx- The former, 
having a wavelength of typically nano meters, will not contribute to the interference of a 
macroscopic length period. Consider now only the latter; its Nx is evidently determined 
by the interaction time, Tw, taken for creating the purely rotational wave train. That is, 
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Tw = N\Ts, Ts = 2-nliOs being the time period of the wave (or the period of its generating 
source), and w,, = 2TTVg/\s. Clearly, for wave interference to occur at a low n number in Eq. 
(46), first of all needs be so large as to be comparable with L. Given m, of the helium 
atom here, is large compared to that of the usual quantum particles such as the electron, an 
(arbitrarily) large Ag (= h/mVs) can be achieved by compensatorily making the flow velocity 
Vg (arbitrarily) small, thus the Tg and L being accordingly large. In experimentation, for 
instance, a large Tg value ranging 30 ^ 600 seconds, with a L ^ 0.0025 cm were achieved."^"^ 
Furthermore, to achieve n = N\ > 2 oi » 2, the Tw should > 2T or >> 2T. 

If the above conditions are met. We can then achieve Lw = > 2L, with Nx > 2. 

In this way, each rotational atomic wave train will wind about each loop at least twice, 
interference will now occur. For the loop satisfying 

I = 2ny, (466) 

where the L will sometimes also write as L2n and n is as before an integer, then, the rejoining 
waves will be in phase: 

= e**'^^i(cos27rn + isin27rn) =${£i). 

This gives an amplified total amplitude: ^{£[) + *(fi) = 24'(£i). Re-writing Eq. (48) in 
terms of circulation, K^ot, we get: 

Krot = <f Vgdi — — i> kgM = — kgL = 



m J m m 

n = 0,1,2,.... (49) 

that is, the circulation is quantized. Or, if L = (2n + l)^Xs, the amplitudes will cancel 
out: ^{ii) + = 0. The amplitudes will partially cancel if L is a fraction of Xg and 

n is small (~ 2), but will cancelled out quite cleanly if n >> 2 in which case the partially 
out-of-phase waves are superposed a large n number of times. The above is a standing 
wave to the observer who rotates with the flow. Since all of the rational atom waves in the 
flow have identical Vg, m, and hence identical A^, and arc also in phase (for they arc set 
to rotation at the same time and they are relatively localized), thus they will all present 
amplified amplitudes at the same circumferences L defined by (46b), and zero elsewhere. In 
the former case, the linearly superposed atomic waves yield a total enhanced fluid density, 
presenting macroscopically discrete streamlines. 

Finally, by Sec. VI, the macroscopic rotation following the vessel above is the only rota- 
tional motion, therefore the total circulation is Ktot = i^rot = nh/m, where n = 0, 1, 2, 

The signal of a measurement, if being in proportion to the total amplitude of the atomic 
waves, will then inform a circulation quantization. 

The physical process described above, which can be termed the circular atomic wave 
self-interference, provides a mechanism for the circulation quantization that has been 
experimentally demonstrated. ^ i-i3,3i 



VIII. THE DYNAMIC STRUCTURE FACTOR ABOUT Qb AS DUE TO THE 
SUPERFLUID BOND EXCITATION 



Corresponding to the total excitation energy s{K) (42), the total dynamic structure factor 
from a (pure) superfluid is 

S{q,cj) = S{q,[£{K) + A{K)]/n) (50) 

r 5ph(g,£(i^)), g<g6-f (a) 

= \ 5ph,6(«, [£{K - Kb) + Ab]/h), \q - gbl < f (b) 
I Sb{q,Ab), q = qb (c) 
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here G represents a narrow vicinity of q^. For the phonon excitation, S'ph(5,w), by Eq. 
(23) the e„'s of the scattering system are discrete. Neglecting phonon-phonon collisions, 
which is to a good approximation as based on neutron scattering measurements (Sec II A) 
and thermal conductance measurements, and as has also been theoretically evaluated by 
Cohen and Feynman,^"^ it follows then effectively 

5ph(g,a;) oc ^ [u){K) - i {e^{K) - e„+i(/f))) , (51) 

where K = q{= — kne\), kne and fc^g being the wavevectors of a neutron before and after 
scattering. 

At Qh = In j a = 6, 6 being reciprocal pscudo lattice vector, by Sec. II B S{q,uj) is due 
to the excitation of a superfluid bond and is thus expressed by (50c). In \q — g^l < cr/2 
where when Sb{q,uj) » Sph,b{^,<l), then (50b) writes Sph.b ~ Sb{q,uj); a switching-over 
region where Sh{q,<-^) and S'ph,6 may have similar weights but would be very narrow, and is 
not to be dealt with here. As Sec. II C reviews, Sb{q,to) as measured by neutrons presents 
a peculiarly high and sharp peak centered at ui = Ag/h = 8.6 °K. Below we formally 
express Sb{q,<-^) through a few combinatory modifications from a crystalline solid scattering 
for which we have the precise description, to elucidate the route via which the superfluid has 
acquired this particular scattering feature, apparently differing from that of either a solid 
or a liquid. 

(a). Let us start with an imaginary superfluid, in which the atoms are arranged in a 
periodic lattice of a simple cubic structure, and the atomic bonding and activation energies, 
denoted by A^"^-^ and A^''^ = A'^'^^ + 5'^''^, are as large as in a solid. For this system, (only 
considering the longitudinal modes) the excitation energy £{q) is regularly given by Eq. 
(26) to be a sinusoidal dispersion, the dotted line in FIG. la, which is zero at qb = b. 
£{(lb) = ^{(Iq) = together with qba = 27r implies that at qb, S'ph(5&,w) = and an elastic 
Bragg scattering occurs. The scattering intensity at a given q is proportional to the static 
structure factor, which (strictly) at zero °K is 

1 " 

[5'ph((7,u;) + Sl^^ {q,uj)\5{uj)d{huj) 

= nSl'^q,Lo)\^=n- (52a) 

The first expression of (52a) is just the Bragg scattering condition; it gives S"^''^(g) = 1 at 
q- (rj — Vj) = qba = n2-K (n = 0, 1, . . . ) and zero elsewhere. The last two expressions of (52) 

generalize S"^{q) to include an energy transfer lo but with cj = for the elastic scattering 
of the present system. The elastic feature is readily understood to result principally from 
the inequality: 

nine ■"s ' (530) 

(5 ~ 100 meV) (1 ~ 10 eV for solids) ^ ' 

where k^c (= irineVne/K), Vac and nine are the velocity and mass of incident neutrons, (b). 
We now reduce A"'' (= Afy - 5v"^) to equal the superfluid bond Ab (= A^ - 6^) ~ 8.6 
"K (0.74 meV), but retain the periodic configuration. Then 



^ ^ne 
2mne 



» A,. (53b) 



Thus At, can readily be excited by neutrons. Then, Eq. (52a) is no longer, in principle, an 
exact description for the present inelastic scattering. However, supplementing the condition 
Eq. (53b) to the conservation relations for the inelastic scattering at q = qb: 
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FIG. 9. Schematic diagram of the relationship between neutron wavevectors fc^ and k'_f before 
and after an inelastic scattering, and the monientuni loss qb that is associated with the energy loss 
As. The Bragg condition is approximately met given that As << (the incident neutron kinetic 
energy) and fc^ — fc^ « As/he^t << k^, therefore k'^. The closed vector- loop ending with 

Qi,^' represents the exact Bragg condition. 

|k;,-k„e|=% = — , (54) 
^(fc;'e-fcne) = A6, (55) 

we obtain that (55) can be adequately approximated by 

^nc — ^ncS (55) 

see FIG. 9. Taken together, (54) and (55)' effectively yield the Bragg condition 

96 ~ 2A;ne sin ^ = — , (56) 
a 

which is equivalent to Eq. (52). In other words, (52a) effectively holds correct for the 
periodically ordered superfluid, but with a; = replaced by a; = = A^/Ti. Modifying 
(52a) accordingly, we have 



= y [5b(g,a;) + S'ph(g,a;)]5(a; - A,in)d{huj) 

= nSl-\q,u^)\^^^jn- {qeqb±cT) (526) 
Re-ordering (52b), we get 

Sb'*'iq,'^)\u=A,,/h = ■^<S''''*(g)U=A,/?i 



e--' ■^'>d{oj- As/n) (q&qb) (57) 



This says that Sl^^{q,Lo) is a modified S^'^^{q), and is maximal when q = qt and ui = Ag/h. 
To re-emphasize, the acquisition of Eq. (52b) is a result of Ab being so negligibly small 
- of the typical energy scale of a phonon compared to the incident neutron energy that 
it causes minimal departure from the Bragg condition; Ab is nevertheless sufficiently large 
to be detectable using the thermal neutrons. This incident is analogous to the thermally 
broadened Bragg scattering - damped by a Debye- Waller factor - in solids. 
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(c). We finally return from (b) to the real superfluid by disordering the atoms. A well- 
defined quantity is now, as with any liquid, the mean interatomic spacing a at which the 
radial distribution function g{r) peaks. As is commonly done with regular liquids^°^~^'^* or 
quasi-crystals,"'"'^^ wc; can rciplacc the site summation of Eq. (57) invovling a (5-function - 
with an integral involving g{r) and meanwhile replace the S{lo — As/h) with h{ui — Ag/h); 
further combining with (50b)-(50c), we get: 



1 

h 



5(>-)->g(r) 



h{uj - As/n), 

qeqb±(T (58) 



here h{co - Ag/h) = {l/y/na) J e"[^('^~^-/'')/'^l'da; sharply peaks around Ag/h and has a 
finite width a. a ought to be basically determined by the fluctuation of the atomic bond 
length duo to the short range ordering. (58) is the scattering fmiction due to the superfluid 
bond excitation. Fixing w ~ ujg, we expect to see S{q,ijj) to peak on the q axis; although, 
there is no experimental data available to us for comparison. Two further remarks can 
be made. First, insofar as the "Bragg" scattering is concerned, it is only meaningful to 
compare the Sb{q,ijj)\u;=u,^ of Eq. (58) with the zero temperature S^'^y{q) of Eq. (52a) and 
not with the finite temperature S'^{q){= S'^ {q,u))duj) of the superfluid; the S"^{q) of 
a solid at zero °K only involves the uj = scattering, whilst the S'^{q) of the superfluid 
contains scatterings of all lu. Second, the correspondence between (52) and (58) not only 
exists in a static sense, but also in a dynamic sense, as atoms in the superfluid, viewed in 
^, are localized as in a solid except that they present no long range ordering. As a basic 
difference from that of a normal liquid, on the other hand, the finite peak width in g(r) of 
the superfluid is primarily due to the disordering of the atomic equilibrium positions but 
not atomic diffusion. 

Through steps (a)-(c), the neutron scattering from the superfluid about qb emerges to 
be an inelastic damped-Bragg scattering due to the excitation of the superfluid bond, 
formally described by Eq. (58). 



IX. THERMODYNAMICS OF HE II 

A. The partition function of the respective fluids 

We consider a bulk of superfluid He II in thermal equilibrium at T, having a volume V 
and number of atoms N. Employing the canonical ensemble of the independent variables 
(T, V, N), the total partition function of the fluid formally is 

Z^ = e^Z^, (59) 

where a = s for the (pure) superfluid and n for the normal fluid of He II, $o is the ground 
state potential energy and e^^o/fesT ^j^g corresponding partition fimction, and Za is the 
total kinetic energy partition function. If the Za is known, the thermodynamic functions, 
the internal energy Ua, free energy Aa, entropy Sa and specific heat Cva can then be 
obtained from the standard relations: 

Ua = ksT = $0 + ksT (60) 

Aa = -kBTlnZa = <^o- ksTlnZa, (61) 
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Sa — ~ 
Cva = 



dT 



d TlnZ, 



dT 



\ dT 



(62) 
(63) 



V,N 



The microscopic theory laid out in the preceding sections enables us to treat the two fluids 
as independent, macroscopic systems, and the practical scheme to obtain in the following 
sections IXB-IXD the respective partition functions and the thermodynamic functions. 



B. The pure superfluid 

Consider a pure superfluid of N helium atoms, which are by Sec. II relatively localized, in 
the normal coordinates can bo represented by N independent, identical and distinguishable 
(see the discussion in Sec. VBc), effectively harmonic oscillators. With the single oscillator 
kinetic energy en{K) of Eq. (23), we get the kinetic partition functions for single oscillator: 



n=0 11=0 
^-hu{K)/2kBT 



I _ g-hu{K)/kBT' 



(64) 



and for A'' oscillators: IlKZph{K). Similarly the partition functions for a single and for 
N2 = 7-/V momentous are as Eq. (64) and ^KM^cph{KM), with Q. in place of w. With the 
above we get 



Zs = ^KZ-p\i{K)J].Km^cp\i{Km)- 



(65) 



Substituting Eq. (64) and its momenton coimterpart in Eq. (65) and in turn, (65) in 
Eqs. (60)-(61), employing the identity \n.\U.KZ-p\,{K)IlKMZci,b{KM)\ = X]^lnzph(i4r) + 
'}2km ^"^ZcphiKM), we thus obtain for the (pure) superfluid: 



f/« = $0 + ksT^ 



E 

. K 



1 dzph _|_ ^ 1 dzcph 



zph dT 



Zcph dT 



Um + I nujV{uj) < n{uj) > + 0„„ (Q) 
Jo 



NkBTT^ 



(60a) 



As = Uso + ksT y" In (1 - e-'^'^/'=«^) V{u;)duj + 0am (f^) 

/ 2/ ln(l 
Jo 



Uso + - 



= UsO + 



el 

NksT* 



(1 - e-y)dy 



yHn{l-e-y)\l 



VD y3 



eu -I 



(61a) 
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and simply with the Ug and Ag above in the last expression of (62), we get 



27^4 / r \ 3 



here U, A and S being each in unit of J. In the above, D{u) = Vlo"^ /2'k'^c\ = 3Nh9'^/kB0^ 
(with three longitudinal modes for each K value) is the Debye density of states of N oscilla- 
tors; 6* = 1^ and 6*13 = ^ T^i^^)^'^ < "('^) >= e'^c./^T.^ is the Planck distribution. 
In the third expression of Eq. (60a) and the third of (61a), y — and yu — j^', for 

yo oo: J^" -^dy = J2^^i V^e'Pydy = X)^! ^ = f^- The momenton terms 0^„, 
and Ogj^ are omitted in the final expressions above. 

= ^o + lNkB{eD+J&D) + Ao, (66) 

8 

is the ground-state internal energy, where $0 the ground-state potential energy; < e„ >o= 
^ksOo, <&Ar>o= ^kB^D and Ag are the zero-point kinetic energies per atom associated 
with the u, ^, and R motions (Sec. II). With the Us above in (63) we get 

with Cv in unit of J/°K. Finally, for the superfluid in a steady flow motion of velocity Vg, 
with the Ss of (62a) we obtain the entropy current 

att'^o ^ t \ ^ 

js = p«v«5,=7VfcB—f^ — v«. (67) 



15 yeoj 

With UsQ = — 7.2°K/atom and ci = 239m/sec as given previously, the above functions 
can be then parameterized. The resulting Ug/NkB, Ag/NkB, and Sg/NkB correspond to 
the graphs in the T < 0.6 °K region in FIG. 10 and these will be discussed in Sec. X. The 

resulting Cy(r) (solid line in FIG. 3a) shows a remarkable Debye behaviour: this has 
been used as a fed-back information in Sec. II A. A least squares fit to the experimental 
data between ~ 0.6 °K using Eq. (74)' gives ci = 241 m/sec and 60 = 20 °K. We 
furthermore observe that, the Debye region in the experimental Cy ends at about 0.6°K, 
which is in rough accordance with the usual T'^ law validity criterion, T < 9d/50 — 0.4 
°K. AnomalousUy, the A profile immediately above 0.6 (1) °K results from the second order 
phase transition and should not be confused with a regular ending of the region. It 
also is anomalous that, at T (=2.17 °K) far below 9d (20 °K), the superfluid has already 
transformed to a normal liquid. 



C. The normal superfluid 

For the normal fluid of A'' identical, distinguishable and weakly interacting classical He 
atoms, the Boltzmann statistical mechanics expressions can be readily written down. The 
total partition function is 

^n = e-^^A (68) 

where the N factory reduces the identical distributions (due to the permutation among the 
N particles) to be counted as one distribution. Uon is similarly defined as Uso of Eq. (66), 
and can in practice be taken as Uno = Ugo, since at zero °K it reduces to the same ground 
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state as for the superfluid. z = YZ^i ,i,=o e"-^*«"*!"*-/''«'^ and = {il + + il)h'^ /SmV^/^ 
are the single atom partition function and kinetic energy. 



v_ 

A3' 



(69) 



where A = h/ ^fhrmkBT is the thermal de Broghe wavelength. Using z of (69) in Eq. 
(68) gives lnZ„ = -Uso/kBT + N\nz - N{\nN - e) = -Uso/keT + N\n{Ve/NA^), where 
e = 2.7183. Placing this in equations (60)-(63) gives: C/„ = Uso + ^NksT (Eq. 60b) and thus 



E = ^NksT, An = Uso + NksTln 



N 



Ve (27i-mfci5T)3 



(Eq. 61b), Sn 



-NkB In 



N 

e5/2 V 



(Eq. 62b), and Cvn = ^Nke (Eq. 63b). 



D. The two fluid coexistent 



1. Two fluid fractions. 

Consider now He II is heated slowly such that at each instant the liquid can be regarded 
as in (quasi) equilibrium; let this start from zero Kelvin where He II is a (pure) superfluid, 
and Vs{0) = V and Ns{0) = N. At a given T, by Sec. II C the fluid is a coexistent of the 
s regions (Sec. II): V,{T) = J2^V,„ 7V,(T) = J^.Ni, and the n regions: Vn{T) = 
Nn{T) = J2i Nni] thus V = Vs{T) + Vn{T) and N = Ns{T) + Nn(T). Furthermore, by Sec. 
II B the s to n conversion involves the conversion of the superfluid bond A;, to the ''van der 
Waals bond" . The probability of the conversion is proportional to the Boltzmann factor, 
l/jgAsn/fesT _ 1] ~ g-Asn/fesT^ bclug as given in Eq. (2a). The density fractions of the 
two fluids of He II then follow to be 

Xn{T) = ^ = x.„e-^ = e-^^'-^'^^\ (70a) 

Xs{T) = ^ = l-Xn = l-e-^^^-^'^"^ (706) 

where Xno = b^bTo jg determined by satisfying the condition x„ = 1 at T = To, To being as 
defined in Sec. IV. Xs and x„ here and the fs and /„ of Sec. IV, are related by: Xg = fs^ 
and Xn = fn^i and are basically equivalent since ^ « ^ w 1. In Eq. (2a), E{T) — ^ksT; 
and put i5vn ~ <5v, then Ab+d^n « 8.6 °K, where Af, = -7.2 "K/atom.^^'^^ Eq. (2a) is thus 
parameterized as: 

A,„ = 7.2 + 1.4+1.5fcBT. (2a)' 

As to To, we here only look at its limits. By the free energy function of Sec. X, the onset 
of n — s conversion would not occur until As < An, which begins at T ~ 3.0 °K (sec the 
context of Eq. 77). On the other hand, To would not be lower than Tx which is 2.17 °K 
from experiment. That is, 2.17 < To < 3.0 °K; the evaluation of Cy in Sec. IX D 2 will 
show, optimally, = 2.6 °K. With To = 2.17, 2.6 and 3.0 °K in equations (70)-(70)', and 
with the Asn in (2)', we obtain three sets of x„(T, Tq) and Xs{T,To) functions (solid lines 
in FIG. 6b). To = 2.6 °K yields indeed the optimum fractions (solid lines 3 and 3' in the 
lower graph) which agree remarkably well with the fs{T) and fn{T) (dashed lines) obtained 
in Sec. IV, and are bounded by the limit fractions with To = 2.17 and 3.0 °K (solid lines 1, 
1' and 2, 2'). 
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FIG. 10. Theoretical thermodynamic functions of the superfluid (s), the normal fluid (n) and 
their coexistent He II (without subscript) (solid lines): (a) internal energies Us, Un, and U as given 
by Eqs. (60a,b), (71)", (b) entropies SsT, SnT, and ST as given by Eqs. (62a,b) and (73)", and (c) 
free energies Aa, An, and A as given by Eqs. (61a,b) and (72)". Circles in (b) represent experimental 
calorimetric entropy data ??. The maximum of A of graph (c) yields the theoretical superfluid phase 
transition temperature T^^'. Ta = 2.17 °K is the experimental phase transition temperature at SVP. 
For the superfluid component of He II, the functions are each given by a canonical ensemble average 
over the states, with a Debye type of density of states, of the independent, indistinguishable single 
phonons, which are produced due to in real space the harmonic oscillations in the normal coordinates 
of the (iV) independent, localized superfluid atoms relative to ^. For the normal fluid, the functions 
are each a Boltzmann statistic mechanics result. The total functions result each as the weighted 
sum of the respective functions of the superfluid and the normal fluid, where the weighting function 
is given by the probability for the excitation of the superfluid bonds. 
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2. Total thermodynamic functions. 



With the Xs{T)^Xn{T) being now known, the total thermodynamic functions of He II can 
be obtained as: 



S — XgSg -\- XfiSji^ 



Cv = XsCvs + XnCvn + -gfi^s - Un), 
'iEntr = Ps'^sSsINUb + PnVnSn/NkB- 



and 



(71) 
(72) 
(73) 

(74) 
(75) 



.1. Between to 0.6 °K. Here equations (70) give Xs(T) ^ 1 and Xn{T) ~ 0, as wc 
have already assumed so when representing the He II in this T region with Eqs. (60a)-(62a) 
and (67); compare FIG. 6b. Hence Eqs. (71)-(75) reduce to Eqs. (60a)-(62a) and (67). 



.2. Between 0.6 (~ 1) to 2.17 °K (T\). Here a drastic conversion from the super- 
fluid to the normal fluid occurs, with both of the fluids in Eqs. (71)-(75) having significant 

weights: < x„ < 1. By Sec. II C the two fluids of He II here aggregate in macroscopi- 
cally large separate regions, thus the lattice statistical thermodynamics results (Sec. IX B) - 
strictly each plus a higher order term Oy{T) {YUg, Ag, ov Ss), and the Boltzmann statistics 
results (Sec. IX C) apply to each fluid regions, across the entire transition T regime except 
in the immediate proximity of the critical point T\ where the atomic dynamics of both the 
fluids become strongly long-range cooperative, for which our present treatment would not 
hold vahd. Using in Eqs. (71)-(74) the thermodynamic functions of Eqs. (60a)-(62a) and 
(67) for the superfluid terms, of (60b)- (62b) for the normal fluid terms, and the Xs and Xn 
of Eqs. (71)-(71)', we obtain: 



U = Ugo-{Ugo-Uon)e 



R T V To ' 



+NkB 



A = UsO- {UsO - Uon)e 



TT^T^ 3 ■ 



159% 



Tin 



Ou{T) 
(71)" 



/ N 



\Ve (27rTOfcBT)3/2 



+Oa{T)}, 



(72)'' 



S = Nkn 



(73)" 




; 1 



kBTj 5 {eDj 



1 



4:kBTj 



-) 



(74)" 

We flrst evaluate Cy. Using in Eq. (74)" the two limiting values 2.17 and 3 °K for Tq as 
discussed in the beginning of this section and A,5„ of Eq. (2b), we obtain two Cy(T) limits 
(solid lines -1 and -2 in FIG. 3b), which fall near to, actually enclosing in the middle, the 
experimental A-profile (circles), therefore defining a region of uncertainty in the theoretical 
Cv{T). Within the Hmits, we find at To = 2.6 °K the Cy(T) (sohd Hne -3) optimally 
matches the experimental Cv data (circles). The remarkably good prediction of Cv and 
similarly Xg {xn) earlier both with To = 2.6 °K, although being fortuitous in principle, can 
however find the following physical explanation. First, an effective To intermediate between 
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2.17 and 3.1 °K is expected to exist as discussed earlier. Second, it can be shown that the 
A-profile is primarily sensitive to the exponent, in which, the only other parameter A^, thus 
Asn, is given consistently by various experiments (cf. Sec. II). 

The first order functions, Eqs (71)"-o(73)" for the coexistent, Eqs. (60a)-(62a) for the 
pure superfluid and Eqs. (60b)- (62b) for normal fluid, can be parameterized, with the values 
for 6d, As, and Tq (2.6 °K) as above and Uqs = —7.2 °K/atom as used several times earlier. 
The resulting hmctions arc graphically shown in FIG. 10, a, b, c, across the whole T range of 
He II; for the entropy there exist experimental data (circles in FIG. 10b), ^^^^^ with which the 
theoretical function S{T) agrees satisfactorily, except for a noticeable suppression about the 
Ta as can be expected. We observe several as-expected features for T < Ty. 1). The internal 
energy and hence the kinetic energy < e„ >= [Ug — Usoj/Nkg, the entropy (at the high T 
end) and the free energy of the pure superfluid are qualitatively lowered compared to those 
of the normal fluid, by < e„ > -Enr = -{Us - Un)/NkB ~ 3.1, {S^ - Sn)/NkB ~ 2, and 
{As-An)/NkB = l°K/atom respectively at T just below 2.17 °K. 2) Prom about 0.6 °K the 
total functions (of the c;oexistent) asymptotically approach to the corresponding functions 
of the pure superfluid with falling T, simultaneously He II approaches to a predominately 
pure superfluid: Xg — > 1. 3). The entropy of the normal fluid from 2.17 K reduces rapidly 
with falling T from Tx, and to a negative function in the regime below ^ 0.6 °K, in which 
it becomes the significant source of the large positive An, or, of the normal fluid phase 
instability. Instead, the localized coherent oscillation scheme of the superfluid produces here 
a relatively high degree of disorder. 4). Us — Uso on the scale of Uso- 5). With Od = 20 °K 
obtained above, we have the zero-point vibrational energy per atom of the u- motion (see 
the discussion below Eq. 66): 

< e„ >o /NkB = 7.5 {"K/atom) (76) 



X. THE SUPERFLUID PHASE TRANSITION AND Ta 

Near Tx, from FIG. 10 we observe that the first order functions of each pure fluid compo- 
nent all suS'er discontinuous changes on transforming from one (normal) component to the 
other (superfluid). However, as the phase transition involves a gradual conversion from one 
(the normal) component to the other (the superfluid), from x„{To) = 1 to a;„(0) = 0, the 
flrst-order total thermodynamic functions of the coexistent. He II, are continuous; this is as 
expected for the second order phase transition here. 

Futhermorc, from FIG. 10 the onset of the n-s conversion can clearly be seen to occur 
near to 3.0 °K, where a lower An switches to a lower Ag: 

Ag{Tx) - An{Tx) < 0. (77) 

With the Ag and A^ given in (61a) and (61b), (77) becomes 

^-19.2{T/enf < (77y 

Since 19.2{T/6d)^ « 1 for the He II temperatures, we can write the left hand side of (77)' 
as q-'^^-HI'/^d) ~ I a,nd re-ordering, we obtain 

rW^r^V^' (78) 
-[\re) ^^^> 

It is noteworthy that this T^^ , as a direct statistical thermodynamic result for the two fluids 
of He II, is basicall}^ identical to the Tx given by London, ^^'"''^ although London's Tx was 
obtained, whilst lacking a direct thermodynamic scheme, by imposing a limitation in the 
thermal wavelength obtained for the normal fluid based on Boltzmann statistical mechanics. 
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Placing now in Eq. (78) the experimental critical density of He II, n = ^ = 2.182 x 10^* 
atoms/m^, and the values for the other constants as previously specified, we obtain 

T^^^ ~ 3.05 "K. (78)' 



However, theoretically the initial n-s conversion under the condition (77), onsetting at T^^\ 
does not lead the total free energy A to maximum until T'-^\ as can be seen in FIG. 10c 
(solid curve). It is the A maximum which resembles the critical instability of the system in 
the second order phase transition, and where a qualitative n — s conversion will occur. This 
critical condition formally writes: 



d^A 

0^2 



= 0, (79) 



which expands as 



- |)(A„ - ^) + - ^ 



+ §f2i^n-A,)=Q. (79)' 

Here the first order term of (A„ — Ag) represents just the condition (77), and the first 
and second derivative terms of [A^ — Ag) are additional; (79)' is a complex second order 
differential equation of An — As- Now instead of seeking an analytical solution of Eq. (79), 
we can readily graphically locate, as we just actually did, the critical temperature from the 
A maximum in FIG. 10c to be: 

^ ~ 2.6 °K. (80) 

Compared to t|^\ T^'^ is seen to be lowered in value in closer agreement with the experi- 
mental value 2.17 °K, although with still a quite large discrepancy. The major error source 
may be the neglect of a strong atomic correlation which anticipates to present already in 
the normal fluid in this T range, in presumably a similar way as in its destination state - 
the superfluid; indeed previous studies^^ have actually shown that the normal liquid He I 
deviates from ideal gas when nearing T\. Correlation will suppress the An value, thereby 
yield a lowered value for Tjj . If assuming across the transition T region A„ is uniformly 
modified to A'^ and hence A to A' and letting the maximum of A' yield the experimentally 
observed Ta, we find A'„ « A„ - 0.7 °K and A! =A- 0.7 °K (dashed curve in FIG. 10). 
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APPENDIX A: DISTINCTION OF CONDITIONS FOR "FLUIDITY" AND 

"DIFFUSIVITY". 

The preceding discussions in effect imply that with the superfluid He II fluidity is not 
accompanied by any discernible atomic diffusion. In combination with the equations of 
motions and their solution in Sees. VB and VC, and, based on the well-understood charac- 
teristics of the liquid and solid systems, we now elucidate that these two kinds of phenomena 
do not need co-occur. Fluidity is a phenomenon that any part of the fluid can move rela- 
tive to another, in a flexible manner and under a small degree of shear stress. In respect 
of motion, the realization of a "flow" only requires the collective translation of all of the 
atoms concerned; there is nothing specific about whether the individual atoms diffuse in the 
same time. In respect of the atomic interaction, fluidity results from the absence of shear 
elasticity, or, the condition that the shear attraction between the atoms is so weak such that 
the system is unable to resist any shearing pressure. This character, in the matrix form of 
the interatomic force F below, is represented by that the shear (tangential) attraction force 
terms - the off diagonal components - are zero (cf. also Eq. (22)): 

Fxx Fxx ~\~ 

F^, 

( liquid at rest); (Al) 

or, 

^ XX ~ XX xy ^ xz 

fp r VP ^ TP TP ^ 

^yx J^yy T JTyy ryz 

p r jp r TP r , jp 

^ zx ^ zy ^ zz ^ ^ z 

(liquid in flow motion.) (A2) 

For the liquid at rest as of (Al), all the off-diagonal force terms are zero, and the diagonal 
components = Fy = F^ consist of the tensile attraction force {F"^) and repulsion force 
i-^aa)^ where a = x,y, z\ the equality of the diagonal terms represents the isotropy of the 
liquid. When in flow motion as of (A2), the shear repulsion force terms, i.e. the off-diagonals, 
are finite and this leads to viscosity. On the other hand, as of Eq. (Al) or (A2), for the 
liquid either at rest or in flow motion, the tensile attraction force terms are finite and this 
represents that the atoms are attracted to each other (true for all liquids). The above are 
to be contrasted to the force condition for an elastic solid: 

TprxTpaTprxTpaTprxTpa 
^ XX I ^ XX ^ xy I ^ xy ^ xz ^ xz 

jp r _i_ TP n TP ^ -L TP ^ TP ^4- TP^ 

■^yx T J^yx^ ^yy "I" ^yy ^yz T -c yz 

TP r , JP a p r , p a p r i p a 

^ zx ^ ^ zx ^ zy ^ ^ zy ^ zz ^ ^ zz 

(elastic body). (A3) 

where, the diagonal terms are generally not zero, whence yielding the rigidity. Diffusion, 
on the other hand, is produced if the gradients of some of the tensile force terms the 
diagonal terms - are finite; e.g. diffusion in x direction results from dF^x/dx ^ which 
leads to a density flux DdF^x/dx; cf. Eq. (13. b). The tensile force terms exist for both solid 
and liquid, this feature corresponding to the characteristically finite values of the diagonals 
in all the above matrices. Hence diffusion can identically occur in both the systems; or, 
they may identically not occur if the kinetic energy of the atoms is low, roughly lower 
than the (tensile) attraction potential barrier between adjacent atoms. It is true that, in 
regular fluids, a trivial shearing elasticity is often accompanied by a high mobility of the 
single atoms. This is nevertheless not so in the superfluid, where the many quantum atom 
correlation (Sec. Ill) produces a contrast in the tensile and shearing interactions, and the 
significantly large tensile attraction prevents the atoms from diffusion. 



39 



^ W H Koosom, Helium, (Elsevier Publishing Comp, Amsterdam, 1942) 

2 W. H. Keesom and J. N. van der Ende, Proc. Roy. Acad. Amsterdam 33, 243 (1930). 

^ J. F. Allen and A. D. Misener, Nature 141, 75 (1938); Proc. Roy. Soc. London A172, 467 

(1939). 

■*P. Kapitza, Nature 141, 74 (1938). 

^ J. D. Reppy and D. Depatie, Phys. Rev. Lett. 12, 187 (1964). 

^ J. B. Mehl and W. Zimmerman, Jr, Phys. Rev. Lett. 14, 815 (1965); Phys. Rev. 167, 214 (1968). 
R. P. Henkel, G. Kukich and J. D. Reppy, in Proceedings of the eleventh International Conference 
on Low Temp Phys, St Andrews, Scotland, Editor J F Allen, D M Finlayson and D M McCall 
(Univ. of St Andrews Press, St Andrews, Scotland, 1968) p 178; BuU. Am. Phys. Soc. 13, 663 
(1968). 

* J. F. Allen and J. Jones, Nature 141, 243 (1938). 
^ F. Simon, Nature 133, 460, 529 (1934). 

^"H. Karncrhngh Onnes, Proc. Roy. Acad. Amsterdam 13, 1903 (1911). 

" W. F. Vinen, Nature 181, 1524-25 (1858); Proc. Roy. Soc. A260, 218 (1961). 

G. W. Rayfield and F. Reif, Phys. Rev. 136, A1194 (1964). 

S. C. Whitmorc and W. Zimnicrniann, Phys. Rev. 166, 181-96 (1968). 

W. H. Keesom and K. Clusius, Proc. Roy. Acad. Amsterdam 35, 307 (1932); Leiden Commun 
219c; W H Keesom and A P Keesom, Proc. Roy. Acad. Amsterdam 35, 736 (1932). 

H. Palevsky, K. Otnes and K. E. Larsson, Phys. Rev. 112, 11 (1958); H. Palevsky, et al., ibid., 
108, 1346-7 (1957). 

J. F. Allen, R. Peierls, and M. Zaki Uddin, Narure 140, 62 (1937). 
^^F. E. Simon and C. A. Swenson, Nature 165, 829 (1950); C. A. Swenson, Phys. Rev. 79, 626 
(1950); 86, 870 (1952); 89, 538 (1953); J. S. Dugdale and F. E. Simon, Proc. Rpy. Soc. London 
A218, 291 (1953). 

^® F. London, Superfluids, Vol 11, Macroscopic Theory of Superfluid Helium, (John Wiley & Sons, 
Inc., New York, London, 1954). 

C. J. Gorter, editor. Progress in low temperature physics, Vols. I- VI, (Interscience Publisher, 
New York, 1955-1970). 

K. R. Atkins, Liquid Helium, (Cambridge Univ. Press, 1959). 

I. M. Khalatnikov, An Introduction to the Theory of Superfluidity, (W. A. Benjamin Inc., New 

York, Amstcrdanm, 1965). 

P. A. Egelstaff, An introduction to the Liquid State, (Academic Press, London and New York, 
1967), Chapter 16. 

W. E. Keller, Helium-3 and Helium-4, (Plenum, New York, 1969). 

J. Wilks, The Properties of Liquid and Solid Helium, (Clarendon Press, Oxford, 1967); An 
Introduction to Liquid Helium, Clarendon Press, Oxford (1970); J. Wilks and D. S. Betts, An 
Introduction to Liquid Helium, (Clarendon Press, Oxford, 1987) 2nd edition. 
2^ A. D. B. Woods and R. A. Cowley, Rep. Prog. Phys. 36, 1135 (1973). 

K. H. Bennemann, and J. B. Ketterson, editors. The physics of liquid and solid helium, Parts I 
& II, (John Wiley & Sons, New York, 1975). 

D. E. Brewer, editor. Progress in Low Temperature Physics, Vol. VII (1978) to Vol XIII (1992); 
D. R. Tilley and J. Tilley, Super fluidity and Super conductivity, 2nd edition, (Adam Hilger Ltd, 
Bristol and Boston, 1986), pp. 8. 

H. R. Clyde and E. C. Svensson, "Solid and Liquid Helium", in Methods of Experimental Physics, 
p. 303-403 (Academic Press, Inc. and Atomic Energy of Canada Limited, 1987). 
P. Nozieres and D. Pines, The Thoery of Quantum Liquids Vol. 2, (Addison- Wesley Publishing 
Co. Inc., 1990). 

R. J. Donnelly, Quantized vortices in helium II, (Cambridge Univ. Press, 1991), Editors: A. M. 
Goldman, P. V. E. McClintock, M. Springford. 
■'^ P V E McClintock, J K Wigmore, Low- Temperature Physics: an introduction for scientists and 
engineers (Blackie, Glasgow and London, 1992). 

A. Griffin, Excitations in a Bose-Condensed Liquid, (Cambridge Univ. Press, Cambridge, 1993). 
^*F. London, Nature (London) 141, 643 (1938). 

L. Landau, J. Phys. (USSR) 5, 71 (1941); 11, 91 (1947). 

R. P. Feynman, Phys. Rev. 91, 1291; 1301 (1953); Phys. Rev. 94, 262 (1954); Prog. Low. Temp. 
Phys., edited by Gorter (North Holland, 1955). 
^''Tisza, Nature 141, 913 (1938 ). 

N. Bogohubov, J. Phys. XI, 23 (1947). 



40 



^'^ M. Cohen and R. P. Feynman, Phys. Rev. 107, 13 (1957); R. P. Feynman and M. Cohen, Phys. 
Rev. 102, 1189 (1956). 

D. Pines, in Quantum Fluids: proceedings of the Sussex University Symposium, Edited by D F 

Brewer (North-Holland, Amsterdam, 1966). 

T. D. Lee and K. Huang and C. N. Yang, Phys. Rev. 106, 1135 (1957). 
H. R. Clyde and A. Griffin, Phys. Rev. Lett. 65, 1454 (1990). 
'^^ P. Debye, Ann. Phys. 39, 789 (1912). 

A concise, clear summary is given e.g. in the review of Ref. 25. 
P. A. Whitlock and R. M. Panoff, Can. J. Phys. 65, 1409 (1987). 

A. Criffin, D. W. Snoke, S. Stringarei, edit. Base Einstein Condesation, (Cambridge University 

Press, 1995). 
*''W. L. McMillan, Phys. Rev. A 138, 442 (1965) 
*® R. Jastrow, Phys Rev 98 1479 (1955). 

S. A. Vitiello, K. Rungo, G. V. Chester, and M. H. Kalos, Phys. Rev. B 42, 228 (1990). 
^° T. McFarland, S. A. Vitiello, L. Reatto, G. V. Chester, and M. H. Kalos, Phys. Rev. 50, 13577 

(1994-11); and Res. therein. 
5^ D. M. Ccpcrlcy, Rev. Mod. Phys. 67 , 279-355 (1995); and rcfs. therein. 

R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals, (McGraw-Hill, New 

York, 1965). 

M. Boninsegni, D. M. Ceperley, J. Low Temp. Phys. 104, 339-57 (1996). 

C. W. Woo, in: The physics of liquid and solid helium, edited by K. H. Bennemann and J. B. 

Ketterson, Parts I, (John Wiley & Sons, Now York, 1975), p. 349. 
®® J. X. Zheng- Johansson, B. Johansson and P-L Johansson, "The superfluidity mechanism of He 

U", submitted to Phys. Rev. (2000); http://arXiv.org/abs/cond-mat/0206339. 
®^ J. X. Zheng- Johansson, et al., private communications with a few specialists in the field (2000, 

2001). 

^'^ J. X. Zheng- Johansson and B. Johansson, in Dynamics in Confinement, edited by B. Frick, R. 

Zorn and H. Biittner, J. de Physique IV 10, Pr7-181 (2000). J. X. Zheng- Johansson, XXH Int. 

Conf. Low Temp. Phys. Abstracts, pp. 296 (Helsinki, 1999. a); Bulletin of the Amer. Phys. Soc. 

- Centennial Meeting 44, 193 (1999. b); Internal report: The Theory of Superfluidity of '^Helium, 

(The Zheng- Johansson Publishing, Nykoping, Sweden, Nov. 15, 1998); Cond-mat/9901125. 

D G Hurst and D G Hcnshaw, Phys. Rev. 100, 994-1002 (1955) 
5^ D. G. Henshaw, Phys. Rev. 119 9-13 and 14-21 (1960). 
™ E. K. Achter and L. Meyer, Phys. Rev. 188 291-300 (1969). 
^1 R. B. Hallock, Phys. Rev. A5, 320-30 (1972). 

W. I. Gordon, C. H. Shaw and J. G. Daunt, J. Phys. Chem. Solids 5, 117-128 (1958). 

W. H. Kecsom and K. W. Taconis, Physica 4, 28, 256 (1937); 5, 270 (1938). 

E. L. Andronikashvih, Zh eksp theo Fiz 16, 780 (1946). 

V. P. Peshkov, Dokl. Akad. Nauk SSSR 45, 365 (1944); J. Phys. Moscow 8, 381 (1944); 10, 389 

(1946); Proc. Conf. Low Temp. Phys. (Physical Society, London, 1947), pp.19. 
W.Y. Tarn and G. Ahlcrs, J. Low Temp. Phys. 66, 173 (1987) 

D. G. Hcnshaw and A. D. B. Woods, Phys. Rev. 121, 1266 (1961). 

A. D. B. Woods, E. C. Svensson, Phys. Rev. Lett. 41 974, (1978); E. C. Svensson, R. Scherm 
and A. D. B. Woods, in Proc. of the Fifteenth Int. Conf. on Low Temp. Phys. (Grenoble-France, 
1978). 

®^ E. C. Svensson, in: Elementary Excitations m Quantum Fluids, pp.59, K. O. Ohbayashi, M. 

Watabe, editors. Proceedings of the Hiroshima Symposium, Hiroshima, Japan, August 17-18, 

1987 (Springer- Verlag, 1989). 
^° W. G. Stirling and H. R. Clyde, Phys. Rev. B41, 4224 (1990). 

O. Plantevin, B. Fak, H. R. Glyde, N. Mulders, J. Bossy, G. Coddens and H. Schober, Phys. 

Rev. B63, 224508 (2001). 

O. Plantevin, B. Fak, and H. R. Glyde, J. de Physiq. IV 10, 177 (2000). 
■^^B. Fak, O. Plantevin, H. R. Glyde, and N. Mulder, Phys. Rev. Lett. 85, 3886 (2000); H. R. 
Clyde, O. Plantevin, B. Fak, G. Coddens, P. S. Danielson and H. Schober, ibid., 84, 2646 
(2000). 

Excitation properties of He II at higher q values than 1.93 1/A have been studied in detail by 
various authors, see e.g. W. M. Snow and P.E. Sokol, J. Low Temp. Phys 101, 881 (1995); T. 
R. Sosonick, W. M. Snow, R. N. Silver and P. E. Sokol, Phys. Rev. B 43, 216 (1991) and Refs. 
therin; B. Fak and K. H. Andersen, Physics Letters A 160, 468 (1991) and Refs. therein. This 



41 



q region involves primarily the translational kinematics of the helium nuclei; this is not covered 
in the discussions of the present paper. 

E. C. Sevnsson, Private communication (1999). 

■^"^ K. R. Atkins and R. A. Stasior, Can. J. Phys. 31, 1156 (1953). 
■^■^K. Skold and K. E. Larsson, Phys. Rev. 161, 102 (1967). 

P. A. Egelstaff, Rep. Prog. Phys. 29 333, (1966). 

J. R. D. Copley and J. M. Rowe, Phys. Rev. Lett. 32, 49 (1974). 
*° J. R. D. Copley and S. W. Lovesey, Rep. Prog. Phys. 38, 461-563 (1975). 

A. D. B. Woods, Phys. Rev. Lett. 14, 355 (1965). 
*2 P. C. Hohenberg, P. C. Martin, Phys. Rev. Lett. 12, 69 (1964). 

As simple cubic, or any other crystal structure, is not an exact representation of the actually 

disordered superfluid, hence the a = 3.3 is only close but not equal to the lattice constant 

Up = 3.6 which one would obtain if assuming a simple cubic structure with the experimental 

density of He II, p = 0.146 g/cm** as input. 

J. Wiebes, C. G. Niels-Hakkenberg and H. C. Kramers, Physica 23, 625 (1957). 
J. D. Wasscher and C. J. Gorter, Physica 18, 329 (1952). 
R. W. Hill and O. V. Lounasmaa, Phil. Mag. 2, 143 (1957). 
*^^D. Greywall, Phys. Rev. B 18, 2127 (1978); 21 1329 (1979) (Erratum). 

F. E. Simon and C. A. Swenson, Nature 165, 829 (1950); C. A. Swenson, Phys. Rev. 79, 626 
(1950). 

E. F. Talbot and H. R. Glyde and W. G. Stirling and E. C. Svensson, Phys. Rev. B 38, 11229 

(1988). 

^° E. C. Svensson, W. Montfrooij and I. M. de Schepper, Phys. Rev. Lett. 77, 4398 (1996). 

J. L. Yarnell, G. P. Arnold, P. J. Bendt and E. C. Kerr, Phys. Rev. 113, 1379 (1959); P. J. 

Bendt, R. D. Cowan and J. L. Yarnell, ibid., 113, 1387 (1959). 
^2 R. J. Donnelly, J. A. Donnelly and R. N. Hills, J. of Low Temp. Phys. 44, 471 (1981). 

G. Carcri, F. Scaramuzzi and J. O. Thomoson, Nuovo Cim. 13, 186 (1959). 

'■^*' F. Rcif and L. Meyer, Phys. Rev. 119, 1164 (1960); L. Meyer and F. Reif, Phys. Rev. 123 727 
(1961). 

T. Ellis, P. V. E. McClintock, Phil. Trans. Roy. Soc. A315, 259 (1985). 
K. R. Atkins and D. V. Osborne, Phil. Mag. 41 1078 (1950). 
'''^F. London, Z. Phys. Chem. (B) 11 222 (1930); Trans. Far. Soc. 33, 8 (1937). 

F. Seitz, The modem theory of solids, (McGraw Hill, New York, London, 1940). 
E. C. Kerr and R. D. Taylor, Ann. Phys. 26, 292 (1964). 

D. G. Henshaw and A.D.B. Woods, Proc. 7th Int. Conf. Low Temp. Phys. editors Graham and 
HoUis Hallett page 539 (North Horland, 1961). 

e.g. E. Merzbacher, Quantum Mechanics 2nd edition, (John Wiley &i Sons, New York, 1970), p. 
107 

See e.g. N W Ashcroft and N D Mermin, Solid State Physics, (Saunders College, Philadelphia, 

1976), pp. 506-8. 

H. Lamb, Hydrodynamics, 6th edition, (Cambridge University Press, New York, 1932). 
H. A. Fairbank and J. Wilks, Proc. Roy. Soc. A 231, 545 (1955). 

Cohen and Feynman's evaluation of phonon-phonon collision time, Phys. Rev. 107, 13 (1957), 
and also that given by Landau and Khalatniko (1949), which corresponds to a peak width in 

the dynamic structure factor < 10~^ °K, should in principle apply to the present phonons, as 
it is chiefly the atomic dynamics scheme for producing the phonons that differs in the present 
theory from Landau and Feynman's theory. 

P. A. Egelstaff, An introduction to the liquid state, (Academic Press, New York, 1967). 

G. E. Bacon, Neutron Diffraction Srd edition, (Clarendon Press, Oxford, 1975). 
G. L. Squires, Introduction to the theory of thermal neutron scattering, (Cambridge Univ. Press, 
Cambridge, 1978). 

C. Janot, Quasicrystals, A Primer, second edition, (Oxford Science Publications, Oxford, 1994). 
The fact that the experimentally observed well defined excitation energy u!{K) makes a com- 
pelling evidence that the distortion here cannot have been caused by the fluctuation in cage 
parameters as we unrealistically proposed earlier. 

C. Kittel, Introduction to Solid State Physics, 6^^ edition, (John Wiley & Sons, Inc., New York, 
1986). 



101 



109 
110 



42 



